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STRONG CYLINDRICAL BENDING OF ELASTIC PLATES 


S. J. Medwadowski,! A. M. ASCE, and K. S. Pister,2 J. M. ASCE 
(Proc. Paper 1692) 


ABSTRACT 


The effect of shear deformation on strong cylindrical bending of an elastic 
plate under uniform transverse load is considered. A system of differential 
equations governing the problem is stated, utilizing the refined theory of 
plates due principally to Reissner and Eringen. The problem is solved and 
graphs are given, from which numerical results are obtainable. The solution 
is compared with the corresponding result of the classical theory. 


NOTATION 


The following nomenclature is used in this paper: 


u, v, W = displacements of the middle surface of the plate 


Wx, Wy = rotations of the normal to the undeformed middle surface 


= Young’s modulus in the x, y plane 


= Airy stress function defined in Eq. (3) 
v = Poisson’s ratio in the x, y plane 

E, = Young’s modulus in the z direction 

. = Shearing modulus in the z direction 


Enh 
~ 12 (1 - v2) 


Note: Discussion open until December 1, 1958. To extend the closing date one month, 
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1. Head, Structural Dept., John A. Carollo, Cons. Engrs., Berkeley, Calif. 

2. Associate Prof. of Civ. Eng., Univ. of California, Berkeley, Calif. 
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stress resultants as defined in Reference (6) 


transverse load on plate surface 


INTRODUCTION 


The classical formulation of the equations governing the flexural behavior 
of thin elastic plates disregards the effect of transverse shear deformation. 
This is characteristic of both the Poisson-Kirchhoff small deflection theory 
and the Féppl-v. Karman large deflection theory. Refinements of the small 
deflection theory to include the above effect may be found in the work of y 
Levy(1) or Reissner, (2) who further considered the effect of transverse 
normal stress. These refinements were also extended to the problem of 
strong bending of sandwich plates by Reissner.(3,4) A recent paper by 
Eringen(5) on the finite flexural motion of visco-elastic plates contains 
Reissner’s static plate equations as special cases, for both small and large 
deflection theories, although no applications to specific problems are pre- 
sented. 

In this paper the refined plate theory is applied to the problem of strong 
cylindrical bending of a uniformly loaded, simply supported strip. The so- 
lution is compared with the results from classical plate theory. The material 
of the plate will be assumed isotropic. 


Fundamental System of Equations 


The fundamental system of partial differential equations governing the 
problem of large deflection of elastic plates including the effect of shear defor- 
mation has been derived and adequately discussed in references cited above. 

In general there are thirteen unknown functions to be determined; introducing 
the Airy stress function F and then choosing F, w, Vx and Vy (see NOTATION) 
as the fundamental unknowns it can be shown that the following fundamental 
system of four partial differential equations must be satisfied: 


F « & ( Way War ) (1a) 


Eb 
tv) G 


P 

1b 
E , z 
W x 7 

= Lu (1d) 

20(1-v) G2 
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where subscripts following a comma denote differentiation and 


N= Nxw,, + Ney Wry + Ny Wy (2) 


The remaining ten unknowns are then given directly in terms of the fundament- 
al four as follows: 


adip) Ge , 


As indicated in the introduction the present theory differs from the classi- 
cal theory in that the effect of shear deformation is considered. Those terms 
which do not appear in the classical theory have been underlined for easier 
comparison. It is seen that one of the equations of the fundamental system, 
Eq. (1a), is identical to the corresponding equation of the classical theory. 


NV, = h V, = h May h (3) a 
re 
4 
ae Dv ) 4a 
AZ 
. 
Sh Gj 
4 
May = + (Yay ™ ao 
q 
Nx v Ny (5a) 
shG,z (5d) 
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The remaining three Eqs. (1b, c, d) contain additional terms. If the effect of 
shear deformation is neglected the terms multiplied by h2 in Eqs. (1b, c, d) 
vanish and the resultant shears Vx and Vy become entirely defined once the 
stress function F(x,y) and weighted displacement w(x,y) become known. 


Cylindrical Bending of an Infinite Strip 


Consider the problem of cylindrical bending of an isotropic infinite strip 
bounded by lines x = 0 and x = a and acted upon by a transverse normal load 
p =const. For a simply supported plate we have the following set of boundary 
conditions: 


ur(o,y) = W(a,y) = 0 
Mx(o, y) Mx (a, y) = (6b) 

either Ww, y) =U 
(6c) 

Mxy(%y) = My (4,7) 
and u = U(4, y) 
ar (0,y) =O 


The solution of the corresponding problem of the classical plate theory can 
be found in reference (6). 


We observe that in this problem none of the quantitites to be found depend 


on y. Also we have My, = 0; Vy = 0; Nxy = 0; hence Nx = const., Ny = constant. 


Then Eqs. (1a, d) are identically satisfied and Eq. (1b), for an isotropic ma- 


E 
terial with G, = av)’ becomes 
D KKK xX = prt N x x x ) Nx LU KKA (7) 


We note next that for the given problem Vx xx = -Nxw,xxx and thus Eq. (1c) 
becomes: 


Vy = 5(i-w) cere (8) 


Thus the order of the fundamental system has been lowered. This is due to 
the fact that we are dealing with a problem of irrotational bending as defined 
by Bolle.(7) To integrate Eq. (7) write it in a different form: 


2 
ur SA. = (9) 


ne 
= 
i 
er. 
| 4 
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2 N » 
(10a) 
tw 
t (10b) 
4 
The general solution of (9) is 
C, X% + C; Stn 2A* + C4 Cos (11) 
DA 
Using boundary conditions (6a) and (6b) we find: 
4 4 
C, = (12a) 
At 
3 
(12c) 
a“ (12d) 
C 


It follows that deflection w(x), maximum deflection w(5) and maximum bending 


moment M,(5) are given respectively by: 


D749 
Lur (x) = fot ~/ x A 
4EDX at* 
(13a) 
c 4 
ee fia) (130) 
> 2 
M, (4%) (13c) 
where 
q 
(a) = | sech | (14a) 
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It is now necessary to find the quantity A. This will be done using boundary 


conditions (6d). On integrating Eq. (5a) and using the first of boundary con- 
ditions (6d) we find: 


—— Ma = dz 
Eh 18) 


On substitution of the proper expression for w,, and performing operations 
indicated we get, after some manipulation: 


< 7 t Z : +9) t aah A (16) 
/¢ ) 4 

which replaces Timoshenko’s formula(8) for classical theory, reference (6). 
This equation can be solved for a given value of the left-hand side (with 
v = 0.3) using appended graph given in Fig. 3. 

The function fj (A) and fg (A) cah be found from graphs given in Fig. 1 and 
Fig. 2 respectively. 

It should be noted that the range of applicability of the figures is limited by 
the requirement that Nx 20. From Eq. (10a) this condition is met if the 
following inequality is satisfied: (for v = 0.3) 


(17) 
where 


re 
/) 

For the case of Nx= 0 (A = 0) the limiting value of f; (A) depends upon the 
span-depth ratio,a@. The value of f; (0) is the shear correction for maximum 
deflection of a simply supported beam under uniform load. 


aa _ be noted that upon setting Nx = -P in Eq. (10a) and by choosing 
- ue —— the Euler buckling load, modified for the effect of shear defor- 


mation, is recovered. 


Example 


To indicate the significance of the results obtained and to illustrate the 
application of the graphs given in Figs. 1, 2 and 3 consider a steel plate strip, 
E = 3 x 107 psi and v = 0.3. Let a = 20", h = 1" and p = 100 lb./in. Then 


Vi9(A) = 2.06 and from Fig. 3 we have A = 0.70. From Figs. 1 and 2, 

f(A) = 0.83 and f2 (A) = 0.70 while the corresponding values from the ciassi- 
cal theory are fj (A) = 0.83 and fg (A) = 0.83. Thus we note that while the ef- 
fect of shear deformation on deflection cannot be observed for the span-depth 
ratio chosen, its effect on maximum bending moment is pronounced. Ex- 


pressions for bending and axial stress are obtainable from Eqs. (13c) and 
(10a) respectively: 


Be: 
: 
~ 
xr 
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0.75 F(x) 


44°) |—— 


( h ha* A (19) 
5 v 


The maximum stress obtained by adding Eqs. (18) and (19) for the case 
considered is 35 ksi as compared to 39 ksi by classical plate theory. A com- 
parison of maximum stresses found by the refined and classical theories for 
span-depth ratios of 20 and 40 is shown below: 


. computed neglecting shear effect. 

It may be noted from the table that maximum stress predicted by classical 
theory is on the safe side. The difference in maximum stress from the classi- 
cal and refined theories decreases as o@ increases, due to the increasing im- 
portance of axial stress and decreasing importance of bending stress. It is 
clear from Eqs. (18) and (19) that the difference occurs in the bending stress, 
while the axial stress is the same by either refined or classical theories. 
While it appears from the results shown that little is to be gained in employ- 
ing the refined theory for normal strength steel plates, its use is justified for 
high strength plates if full advantage of high strength is to be taken. 


CONCLUDING REMARKS 


The problem of strong cylindrical bending of an elastic plate-strip under 
uniformly distributed load has been solved in this paper. The behavior of the 
plate was considered in light of the refined plate theory of Reissner, which 
includes the effect of shear deformation. Graphs are given to simplify 


bis ASCE 1692-7 
= 
4 bend. (18) 
“a 
$$$ 
| 
| 2 = 20 x = 40 
60 23 26 52 5 a 
100 35 39 77 83 
200 57 66 110 120 
3 
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numerical solution of problems and results are compared with those of the 
classical theory as given by Timoshenko. It is shown that, as might be expect- 
ed, the effect of shear deformation is more pronounced for plates whose thick- 
ness is large compared to their span. In the case of maximum bending 
moment, however, the effect of shear deformation appears to be large even 

for fairly thin plates. The importance of the refined theory for efficient de- 
sign of high strength plates is indicated. Plate strips with other types of 
boundary conditions may be treated by the method described. 
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DYNAMIC ELASTO-PLASTIC RESPONSE OF RIGID FRAMES 


Frank L. DiMaggio,! A. M., ASCE 
(Proc. Paper 1693) 


INTRODUCTION 


The procedure for determining the elastic response of structures to 
transient loads in terms of the elastic normal modes is well known.(1) As- 
suming ideal plasticity, Bleich and Salvadori (2) recently proposed a method 
for obtaining the response beyond the elastic range in terms of elasto-plastic 
normal modes. They applied the method to free-free beams acted on by im- 
pulsive loads. 

In this paper, the method of Bleich and Salvadori is applied in general to 
two hinged and fixed rigid frames acted upon by a time dependent concentrated 
force at the knee and in particular to a two hinged frame (for which the elastic 
normal modes are known(3)) and a linearly decaying forcing function. 


General Analysis 


The structure and loading considered in this paper are shown in Fig. 1, 
where the ends A and D may be either both pinned or both fixed. 


First Elastic Stage 
The forced elastic deflections of the columns and beam (see Fig. 2) may be 
represented by 
Uy 
~ (X ,t) = ) t 


Note: Discussion open until December 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE, Paper 1693 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 3, July, 1958, 
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Fig. 1. Loaded Rigid Frame 


(x",t) = ay(t) oy(x") 
isl 


in which $; are the antisymmetrical normal modes of the frame (the sym- 
metrical modes are not excited by the load considered) and qj, the generalized 
coordinates, satisfy the differential equations 


2 
+ wy a = (2) 


Fig. 2. First Antisymmetrical Mode 
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where w; are the natural frequencies, Cj the generalized masses: 


and Q; the generalized forces due to the applied force F(t): 


= F(t) (4) 


In Eq. (2), dots denote differentiation with respect to time t. If the structure 
is initially at rest 


a,(0) = 4,(0) = 0 


and the general solutions of (2) are 


t 


a, (t) = sin w,(t - t) dr 


When, at time t,, the elastic moment at the knees, computed from 


M, (Lt) = - a(t.) (LY) 


1=1 


where t, denotes the time variable during the elastic phase reaches the ca- 
pacity value Mo, plastic hinges are assumed to form there and the first plastic 


stage begins. In Eq. (7) and hereafter, primes indicate differentiation with 
respect to x. 


First Plastic Stage 


On the assumption of an ideal plastic hinge, the moment at the knees re- 
mains constant and equal to My throughout the elasto-plastic stage. If the 
capacity moment is not reached at any other section during this stage, the 
motion of the frame may be determined by considering two separate elastic 
systems as shown in Fig. 3. 


Consider the translational motion of the beam. Newton’s second law yields 


s(t) = re) - (8) 


where S(t) is the shear at the column top as shown in Fig. 3. Thus the effect 
of the mass of the beam on the lateral motion of the column may be reproduced 


by adding a concentrated mass .5UVLY to the column at x¥ = LY, This permits 


the use of the systems shown in Fig. 4 to determine the elasto-plastic 
response. 


3-3 
= & dx +4 (o,)° dx je (L) (3) 
3 4 
Li. 
(6) 
(7) 
q 
por “a 
3. 
fed 
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— 
F(t), B 
M, 


Fig. 3. Column and Beam Elastic Systems 


Denoting the modes of free vibration of the column with concentrated mass 
(Fig. 4a) by Y , &), the lateral deflections may be expressed as* 


Fig. 4. Column and Beam Elasto-Plastic Systems 


| 
f *The symbol ) will be used to denote » for the two hinged frame, where i = 0 


i i=0 
corresponds to the rigid body mode, and » for the fixed frame. 
i=] 
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Yp (x = V(x") 
i 


where the generalized coordinates PP satisfy the differential equations 4 | 
u a 
eu u,2 Py 
Py + py = = (10) = 
In Eq. (10), 2}! are the natural frequencies associated with the ¥}, a are the pO 
generalized masses: 
i (11) 
and P} are the generalized forces associated with the external force and a 
couple: 
u 
P t ss u 
Similarly, the lateral deflection of the simple beam of Fig. 4b is given by 4 
Vv, _Vv Vv 
t) =) \ 
i=1 
where VY are the antisymmetrical modes and the generalized coordinates 
p} (t) satisfy the equations 
5 pv 
i 
Pp + (25) pes (14) 
| 


ay = (ets (15) 


P,(t) = mM, vy (0) (16) 


Initial conditions on the elasto-plastic generalized coordinates p? and py 
are obtained by requiring continuity of displacements and velocities at the 
instant t, when the transition from the elastic to the elasto-plastic stage oc- 
curs. 

Considering the column first, 


yt.) = yp(xyt,) (17a) 


u eu u 
Yp(x rt) s Yp(x st) (17b) 
Using Eqs. (1a) and (9), Eq. (17a) becomes 


i=] i 
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Eq. (18) holds at x¥ = LY: 


(19) 


Multiplying both sides of (18) by puypp (x") and integrating along the column 
yields 


2 
i=l 
Both sides of Eq. (19) are now multiplied by UY LY vu (LY): 
i=l 
Adding Eqs. (20) and (21) results in 
L 
Y (x) V(x) dx + YL" (L") 
4 (22) 


L 
u u u lv 
Making use of Eq. (11) and the orthogonality condition for the modes /}, 


50x, 
u 


L 


Eq. (22) reduces to 
u u u u 1 
In the same way, Eq. (17b) leads to the condition 4 
1 
i=1 
Considering the beam, the requirement that : 


st.) = yp(x",t,) (27a) 


= I(x", t,) (27b) 


and the orthogonality condition for the modes oy : 


he 
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u 


L 

Vy (x) dx = 0 (28) 

leads to the relations 
Py(ty) ) a, a,(t,) 
= 1 
% 
— in which a! is given by Eq. (15) and 
v 
L 


Let 


t-t ew f 
° 


s Then the solution of Eqs. (10) and (14) for i 4 0 may be written as 4 


M 
oi u 
- ———_ / sin -7T) ar 


u 
¥y(L") Tt 


u 
+t.) sin - 7) ar 
I (0) 
P(r) = cos air + sin Q’r 4 
Vv i 

The initial values (i.e., values at t = tg or T = 0) appearing in Eqs. (33) and 4 
ahi (34) are computed from Eqs. (24), (26), (29) and (30). Since : 


Tt 


0 = (2 cos ay 


Eqs. (33) and (34) simplify to 


are 
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B, (0) 
Py(t) = py(0) cos T+ sin t 
4 
- (1 - cos a'r) 
(22) 
Wit") 
aa? [re +t,)ein (35) 
B,(0 
Pi ( t) = p,(0) cos + sin 
aM, (0) 
* (2 - cos 25 t) (36) 
ay (24) 


In the case of a two hinged frame, a rigid body rotation of the column is 


possible in the plastic stage. Denoting this mode of motion by yu , it may be 
expressed as 


Vo" 


(37) 


The differential equation governing the generalized coordinate p¥ associated 
with this mode is obtained by setting 2U = 0 in Eq. (10). Eqs. (11) and (12) 
remain valid. Thus, making use of (37): 


M 
_ F(t) 
p= 38 
° ( ) 


The solution of (38) is 


u 7 F(T + t,) Mo aT 
t) = + if ot" (39) 


Since 


Eq. (39) may be written as 


vA 
M 
u u ° 2 
p(t) = p (0) + p(O)r - a dz 


Except for Eq. (35), all other equations derived for the p? are valid for the 
rigid body mode. 

During the first elastic stage, the beam remains perpendicular to the 
columns at the knees of the frame. After a plastic hinge forms, this constraint 
no longer exists and energy is dissipated as work is done by the couples with 


ane 
{ 
Tt z 
M, M, 2 
— 
u u 
L aL 
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moment M, in moving through a plastic angle. The plastic angle Op (see Fig. 
5) is given by 


= (LY) - (0) 


(41) 
which becomes, using Eqs. (9) and (13), 


= ) - ) vy (0) (42) 
i i=l 
Plasticity at any section is characterized by the condition that the rate of 
plastic work is positive, i.e., 


Wi, = > 0 (43) 


Assuming that the capacity moment M, has not been exceeded at any other 
section of the frame, the time T, at which the first elasto-plastic stage ends 
and the second elastic stage begins is determined from 


6,( t,) = O (44) 

It is to be noted that the bending moment at the knee of the frame may not 
be computed from an expression like (7) since each ¥¥"(L¥) = 0. For points 
near the knee, the convergence will be very slow. In ‘order to find bending 
moments near the top of the column, therefore, the equilibrium of the column 
under the action of external and inertia forces is considered. This will be 
illustrated here only for the case of a two hinged frame. Fig. 6 is a free body 
diagram of the column, excluding longitudinal forces, whose effect on bending 
is neglected in this study. 

Equilibrium of forces in the y direction requires that 


H= 3 - - ax 


Fig. 5. Plastic Angle 
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Fig. 6. Free Body Diagram for Column 


The bending moment at any section x", found by considering forces acting to 
the left of the section, is 


u 
x 
u u 
M(x) = Hx + 2° dt (46) 
which becomes, using (45), 
u 
L 


1 .u,,u sadly 
u 


x 


Introducing the series expansion for ¥p into (47), i.e., 


gx") =) (48) 
and letting 
ax = (49) 
x" 
x Wy (x) ax = 


Eq. (47) becomes 


i 
+ y(t") (51) 


At the knee (x" = L¥) it is easy to show that (51) reduces to Mp(L¥) = -Mo. 


1693-10 EM 3 July, 1958 “ 
A o 4," 
ty" 
a 
: 


ELASTO-PLASTIC RESPONSE 


Second Elastic Stage 


In the second elastic stage, the displacements are again given by Eqs. (1) 
with new generalized coordinates q; satisfying Eq. (2) and initial conditions 
determined from 


U 
»0) 


yp(x", ; 0) = (52) 


40) = yp(x", = 5) (53) 


in which a new time variable has been introduced in the second elastic stage: 
t* = t - (t, + (54) 
Consider the first of Eqs. (52) and (53) and the particular case of Eq. (52) for 


xU = LU, Substituting the appropriate series expansions, these become 


i=l q 
= 
q 
i 


The following operations are now performed: (1) both sides of Eq. (55) are 
multiplied by “4@"(x;,) and integrated along the column,(2) both sides of Eq. 


(55a) are multiplied by SeYLYOU (LY), (3) both sides of Eq. (56) are multiplied 


by “YM (xY) and integrated along the beam. 
If the three equations resulting from these operations are added up using 
the orthogonality condition for the elastic modes, i.e., 
u 


L 


the result is 


k k=l 


where Cj, ik and j, are given by Eqs. (3), (25) and (31) respectively. Sim- 
ilarly, using Eq. (53), initial conditions on the §; are obtained as 


2 ui 
(0) = oy +) (59) 
k k=1 
Using the initial conditions obtained trom Eqs. (58) and (59), the solution of 
Eq. (2) is 
(0) 


a,(t*) = a,(0) cos w,t* + 


sin w,t* 
i 


o sin a, (t* - T) ar 


ig 
q 
= 
tai 
wee 
ber 
4 

a 
a 
f 
ve 
t# 
(60) 
hy 
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Eq. (7) can no longer be used to find the bending moment at the knee be- 
cause the series obtained by differentiating the series expression for yp twice 
with respect to x does not converge to yf}. at the knee. This can be seen by 
realizing that the series obtained by differentiating each term once does not 
converge to yt, since this series of derivatives can never lead to a plastic 
angle at the knee. At points near the knee, the series of second derivatives 
will converge slowly to the curvature, 

Using inertia forces, a series expression similar to Eq. (51) is obtained. 
For the numerical examples considered in section III, however, this series 
did not seem to converge either. 

It therefore appears that it is not possible to compute the moments near the 


knee after the first plastic stage and thus impossible to know whether a plastic 
hinge reappears there. 


Two Hinged Frame and Linearly Decaying Load 


In this section, the results of section II will be particularized for the case 
of a frame supported by smooth hinges and a linearly decaying load, i.e., 


F(t) 


t 
t<t 


e 


u 


, t>t, (61) 


First Elastic Stage 


In a recent study, DeHart(3) obtained the modes and frequencies for a two 
hinged frame and determined the response to the loading given by Eq. (61) up 
to the end of the first elastic stage. His results will be summarized in this 


subsection. The elastic frequencies wj associated with the antisymmetrical 
modes were obtained from the equation 


k,L 
cosh + - cos —5— sinh (64) 


ny = sin 


J 
‘ 
o 
4 4 
where 4 
o (ky) (uy) 
i} 
and it 
viv 
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Vv i 
Ns = 2 sin 5 sinh 5 (65) 


ny = sin cosh - cos kL" sinh 


(66) 


u u u 
No = 2 cos cosh (67) 


u u u 
13 = 2 sin sinh (68) 


The antisymmetrical modes are given by 


in which 


u 
= sin + D, sinh 


v 


i 


i 


(cos - cosh + sin k 


u 
sin + D, sinh 


(72) 


Using the modes given by Eqs. (70) and (71), the generalized mass Cj of 
Eq. (3) is obtained as 


| 
a 

Fy 

a 

= 

k,L 
v i 

i i 2 
bss 

tan k,L 7 
i a 
2 vv 

‘ (70) 
2 

1 

uu k,L / v 

- cos k + — sin k 
2 uv u i 4 

A 

4 

4 

1 | tan k,L 
/ = | | 
hy 

| 

ths 

| 
& 4 
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2 
ein + D, sinh 


1 u u 


(cosh gin KUL" - sinh cos KUL") 
u i i i 


(73) 


and, for the forcing function given by Eq. (61), the elastic generalized coordi- 
nates of Eq. (6) become, for t < te, 


i 


t | (74) 


First Plastic Stage 


The normal modes of the column with concentrated mass (Fig. 4a) satisfy 
the differential equation 


u od u u 2 u 
and the boundary conditions 


(75) 


= = = 0 


2 

(i) + =0 (76) 
In the usual manner, it is found that solutions to Eq. (75), satisfying the 


conditions (76), exist only for discrete values of the frequency determined 
from 
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ky 
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+ 
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k,L 
i 
1 1 i 
+ ——— |= /1 - --———— _ | + — sin —— cos —— 
= te \ v,Vv 2 2 
[ \ k,L k,L 
sinh 
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u u 
(cot - coth = rar (77) 


By 


u u 
4 
78 
i EIU (78) 
Corresponding to the frequencies Q} found from Eq. (77) (which has a solution 
2, = 0) the normal modes are 


(79) 


(80) 


For the simple beam of Fig. 4b, the well known antisymmetrical modes and 


frequencies are 
Vv ax” 
= sin i — 
i LY 


2 
Vv 
2 
Using the expressions (69) and (70) for the elastic modes and (77) and (82) 


for the elasto-plastic modes, the generalized masses of Eq. (11) and (15) and 
the coupling constants of Eqs. (25) and (31) are obtained by integration as 


u 
u 


a, 2 (83) 


2 
yon 2 u u 
+ 3 sin® B,L coth ByL 
sinh 


ay = - 


Ps 


uu 
a = (in a + D, a) + [sina - cos + D, 


(a cosh a - sinh a) | 


ASCE 
where 
i 
it 
- 
y. = sin 8.x + ——— sinh , 1>0 
i i i 
sinh 6, 
: 
¥ 
\ - 
F 4 
2 
‘ 
q 
if 
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= gin b (sin a + D, sinh b) 


3 
+ at 2a” sin b (D, cosh a - cos a) 


+ a> (cos b + sin b coth b)(sin a + D, sinh a) 


+b? (cos b - sin b coth b)(sin a - D, sinh of j>o (87) 


where 


J (88) 


(89) 


For the forcing function of Eq. (61), the integrals in Eqs. (35) and (40) be- 
come, for t < tg, 


F 
+ t,) sin -T)@e= - - cos ait) 
ive 


sin ta 
(90) 


u 
F t 
1 0,2 


If bending moments near the top of the column are desired, the integrals of 
Eqs. (49) and (50) must be evaluated. From Eqs. (79) and (80) these are ob- 
tained as 


(92) 


(93) 


u (94) 


at 
k,L 
Vv mL Vv \ 2 
a 
( kL 2 
| bi 

ha 
| 
1 
u x" 
u 
Ks (x )= = | - cos + (cosh -1)| ; 
sinh 
i>o 
3 
P, u x" 
= 

3L 4 
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P,_u 1 uu uu uu 
)= sin Bx ~ B,x cos 


u 
(6;) 
sin 95) 
einh 
Second Elastic Stage 
: The integral in Eq. (53) becomes, for the force function of Eq. (61) 
Fo lee (a os + sin - (96) 
@,t 1 @ 4 
Fy 2. 
E (1 - cos @, += ein |; t>t, (97) 
ive 1 
where 
(t, + (98) 
Numerical Examples 
Using a forcing function given by Eq. (61) with time constant 
= .5 sec 
and two values of the maximum force 
(1) F, = 345600 lbs. (2) F, = 57600 lbs. 2 


response calculations were made for a two hinged steel frame with the follow- 
ing properties: 
LY = 14s in. LY = 360 in. 


“ 


1" = 39% in. 1” = 1043 in. 


M, = 2.3 x 10° in. lbs. q 
This is the same structure used by DeHart, (3) who calculated the response in 
the first elastic stage to the higher of the two loads considered here. 

ee Throughout the calculation, infinite series were approximated by the sum 
ae of the first five terms. All numerical work was done with an ordinary desk 


computer. 


Elastic Modes and Frequencies 


The first five frequencies were obtained from Eqs. (62) and (63) and 
substituted into Eqs. (69)-(73). Useful results are summarized in Table 1. 


“dies 
] 
q 
| 
| 
| 
| 
q 
g 
q 
a 
al 
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Elasto-Plastic Modes and Frequencies 


The column and beams frequencies were obtained from Eqs. (77), (78) and 
(81) and substituted into Eqs. (79), (80), (83)-(89). These results are shown 
in Tables 2, 3 and 4. 


Response to Load with F,, = 345600 lbs. 


First Elastic Stage.—The deflection and bending moment at the knee, ob- 
tained by substituting the expressions (69) and (74) into Eqs. (la) and (7), were 
calculated as functions of time using the constants of Table 1. The time at 
which the bending moment reached the capacity value My was found to be 
te = .044 seconds. In Figs. 7 and 8 these results are plotted. 

While it was found that all modes other than the first contributed almost 
nothing to the deflection throughout this entire time range, the third mode was 
at least as important as the first mode in determining bending moment for 
early times. This is illustrated in Fig. 8. 

It should be noticed that the deflection at the top of the column when the ca- 
pacity moment is reached is almost exactly the same as that which would oc- 
cur if the load were applied statically. This is, however, not generally true. 
For a different frame or forcing function, the static and dynamic deflections 
at incipient plasticity might be quite different. 

Table 5 shows the values of the generalized coordinates and their deriva- 
tives at the end of the first elastic stage. 

First Elasto-Plastic Stage.—Eqs. (24), (26), (29) and (30), when used in 
conjunction with Tables 2, 3, 4, and 5, furnish the initial values of the plastic 
generalized coordinates and their derivatives. These are recorded in Table 6. 

It should be pointed out that the series expressions for p}'(t,) and p}(t,) 
do not converge well enough to justify the number of significant figures pre- 
sented in Table 6. For the higher values of i the first figure is in doubt. The 
numbers presented in these cases are merely the result of adding the five 
terms of the series, each of which is accurate to the number of decimal places 
in the sum. The spurious nature of these numbers does not, however, effect 
any further results appreciably. 

It is interesting to see how well Eqs. (17a) and (17b) are satisfied at the top 
of the column using the values of the generalized coordinates of Tables 5 and 
6. In conjunction with Tables 1 and 2, it is found that 


yu(L",t,) = 1.970 + .O14 + .000 + .000 +.000 = 1.964" 

yp(L' st.) = 2.028 - .039 - .003 - .001 + .000 = 1.985" 

J_(L yt.) = 86.69 - .82 + .56 + .20 - .05 = 86.58"/sec 

j(b'st,) = 89.41 - 2.61 - .09 + .11 - .02 = 86.80"/sec 
Also of interest is a comparison of slopes at the knee: 


' ' 
ye (h'4t,) = yy, (0,t,) = .004266 + .001705 - .000002 - .000003 - .000002 


= .005964 rad 
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Table 3. Values for ay 


- 
-OL446 
-5651 

-2.036 


2.8225 
-88.122 
32.501 
3.3601 


Table 5. Generalized Coordinates and Derivatives 
at the End of the First Elastic Stage 


‘ 
x 48.59 1.91 .248 .072 | 
2 2h -.1309 .0369 
3 -3.964 2h .06 5.400 -1.374 3 
4 1.326 -4.965 | 22.68 6.478 ae 
5 -.6614 2.398 | -5.177 22.08 7.095 2 
Table 4. Values for a” » 
ij 
1 34904 -10336 043600 .022323 
2 -1891 -05312 -02222 -011348 
3 -7.6812 | -1.9711 | -.6132h | -.41388 
-36.51 9.0360 2.9125 1.4009 
5 7.818 40.16 -9 .6864 -3 4695 
( 3.172 | -.06724 -001715 -003371 | -.0003172 
a | 139.6 3.845 -2.426 1.386 - 4480 
| 
: 
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Table 6. Initial Values of the Plastic Generalized 
Coordinates and Their Derivatives 


Table 7. Elasto-Plastic Generalized Coordinates 
at the End of the First Plastic Stage 


Table 8. Generalized Coordinates and Their First 
Derivatives at the Beginning of the 
Second Elastic Stage 


1693-21 
| 2.028 89.41 2 | .2k85 -.529 
1 179 12.0 4 -01593 3498 
2 | -.0284 | -.79 6 | .005645 .8322 
3 0099 | -1.344 || 8 | .002k72 -1306 
4 10 .001252 07585 “ie 
2 o | 35.16 12.0 2 .0730 3.92 . 
1 O42 5 4 02062 1.30 
2 .0309 -.84 6 .006161 
3 0094 8 002493 113 F 
a -.0051 1. 10 .001273 -O48 
| 
4,(0) | 54.60 | -3.980 1.22 | -.5691 q 
a,(0) | 19-2 -3.7 2 -2.3 .62 
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Deflection at top of Colum, inches 


First Elasto-Plastic 
Stage 


= 345,600 lb 


) Time t, seconds 
| 
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| | 


Figure 7. Relation Between Deflection and Time at Top 
of Column (F, = 345,600 lbs. ) 
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(1",¢ ) = .01408 - .00446 - .00132 - .00068 - 00044 = .00718 rad 
Pp ° 


(Ost) .004337 + .000556 + .000296 + .000175 + .O00109 


-005471 rad 


It is seen that convergence in the elasto-plastic case is slow. As the 
plastic stage proceeds, however, this convergence improves as the rigid body 
mode becomes more and more important. 

If Eqs. (90) and (91) are substituted into Eqs. (33) and (34), the plastic angle 
may be determined as a function of time from Eq. (42) by substituting the ap- 
propriate quantities from Tables 2, 3, 4 and 6. The result is plotted in Fig. 9. 
Except for the small values of T, the contribution of the rigid body mode is by 
far the most important. The time 7, when the first plastic stage ends is de- 
termined analytically from Eq. (44). From Fig. 9, it may be obtained graphi- 
cally as T, = .240 seconds. The maximum plastic angle is seen to be .238 
radians. 

From Eq. (9), the deflection at the top of the column is computed. This is 
plotted in Fig. 7. Again, the major portion of the response is from the rigid 
body mode. 

At time To, the end of the first plastic stage, the elasto-plastic generalized 
coordinates are given in Table 7. 

A check was made of the bending moment at the center of the column during 
this stage using both Eq. (51) and one analogous to Eq. (7). The results, con- 
sidering the convergence problem, were in fair agreement and established 
clearly that the bending moment at this section was well below the plastic ca- 
pacity value throughout the stage. Thus, no plastic hinge other than the one at 
the knee is generated, in accordance with the assumptions made. 

Second Elastic Stage.—Eqs. (51) and (52) together with quantities in Tables 
1, 3, 4 and 7 are used to determine initial values of the new elastic generalized 
coordinates and their first derivatives. These are listed in Table 8. 

Substituting Eqs. (96) and (97), as appropriate, into Eq. (53), with the neces- 
sary values from Tables 1 and 8, Eq. (1) may be used to find the deflection at 
the top of the column as a function of time. This is plotted in Fig. 7. Note 
that the assumption has been made that plastic hinges do not reappear. Again, 


the first mode was by far the most important in the series for the deflection 
at the top of the column. 


Response to Load with F,, = 57600 lbs. 


Figs. 10 and 11 and 12 show the results of this calculation. 

The plot of Eq. (42) for T > 0 (not considering condition of Eq. (44)) shows 
the influence of the higher modes. Since the curve has two stationary points 
before it begins to decrease monotonically, two elastic and elasto-plastic 
stages might have to be examined to get the maximum plastic angle. It should 
be noted, however, that the first dip in the curve is of the same magnitude as 
the truncation error apparent at T = 0, so that the first stationary point may 
be unreal. Since, in any case, it appears impossible to determine when elastic 
stages other than the first come to an end, it seems appropriate to approxi- 
mate the response in the first plastic stage by the rigid body mode only, giving 
the dotted curve shown on Fig. 12. Thus the stage ends, approximately, at 
To = .160 sec. and the maximum plastic angle is .0252 radians. 
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Figure 10. Relation Between Moment and Time 


at Top of Column (F, = 57,600 lbs.) 
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First Five Modes 
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Figure 11. Relation Between Moment and Time at 
Top of Column (F, = 57,600 lbs.) 
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Five modes 
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Figure 12. Relation Between Plastic Angle and Time 
at Top of Column (F, = 57,600 lbs. ) 
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ELASTO-PLASTIC RESPONSE 
CONC LUSION 


It appears that, for engineering purposes, a good approximation of the maxi- 
mum plastic angle which develops at the knee of a two hinged frame under a 
concentrated dynamic load applied at the knee may be obtained by considering 
the first elastic mode and the first (rigid body) elasto-plastic mode only. 
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SYNOPSIS 


Plastic collapse of a partial circle cylindrical shell simply supported at 
each end with the longitudinal edges free is investigated for the case of uni- 
form radial loading. The principles of Limit Analysis applicable to a plastic- 
rigid structure are utilized to develop upper and lower bounds for the collapse 
pressure parameter, P. Bounds are obtained for a wide range in shell shape 
parameters with the deviation between the mean of the bounds and either 
bound varying between 2 and 25 per cent. 


NOTATION 


The letter symbols are defined where they first appear and at the end of 
the paper. 


INTRODUCTION 


In structural engineering, recent years have seen a more frequent adoption 
of cylindrical shells spanning between transverse supports in order to ac- 
complish the roofing of large unobstructed areas. Such a typical shell roof is 
shown in Fig. 1. Current practice provides for the design of such shells on 
the basis of the elastic theory as developed by Fligge, Donnell, and 
Timoshenko.(1,2,3) The practical design of shell roofs by this elastic theory 
has in turn been greatly simplified by the tabulation in ASCE Manual No. 31(4) 
of values of the controlling stress resultants and deflections for various load- 
ing conditions and shape parameters. By this method, the design of cylindrical 


Note: Discussion open until December 1, 1958. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1706 is 
part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. EM 3, July, 1958. 


1. Chief of Design Branch, U. S. Army Engr. District, New York, N. Y. 
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shell roofs using working stresses less than the yield point can be readily ac- 
complished. 

However, for information as to the ultimate load carrying capacity of such 
structures, the theory of plastic collapse must be employed. Limit and plastic 
analyses of full circle cylindrical shells under axially symmetric loading have 
been accomplished by Drucker, Hodge, and Onat.(5,6,7) For shells of revo- 
lution under axially symmetric loading, Onat and Prager have used limit 
analysis to determine the collapse loading. 8) However, collapse loads for 
partial circle cylindrical shells with resultant non-axially symmetric loading 
are not available, and this paper concerns itself with the development of such 
results. 

Several important distinctions between the plastic collapse of full circle 
shells and partial circle shells are to be noted. First, for collapse of the full 
circle shell to occur, the entire shell must be plastic. Although precise col- 
lapse modes are not derived herein, displacement patterns which are utilized 
to develop effective upper bounds on the collapse pressure tend to indicate 
that collapse occurs without the shell being fully plastic. Secondly, the sym- 
metry about the axis of the cylinder which prevails for the full circle no long- 
er is applicable to the partial circle with regard to the resultant stresses and 
deflections. Lastly, four controlling stress resultants must be considered for 
the partial circle, namely, three membrane stresses and the transverse 
moment whereas three stress resultants are predominant in the case of the 
full circle shell. 

In this paper, a linear approximation to the Mises and Tresca conditions of 
plastic yielding is used and the material of the shell is assumed to be a rigid- 
perfectly plastic material. Thus, of the total strain rates the elastic com- 


ponents are neglected and only the plastic strain rates are considered signifi- 
cant. Likewise, by this assumption the effect of strain hardening is neglected. 

A precise solution for the collapse of a rigid-plastic partial circle shell 
would define the zones of plastic action, give the collapse load and the dis- 
placement rate pattern and meet the following requirements. 


a. Throughout the shell all equilibrium and stress boundary conditions 

must be satisfied. 

Within the plastic zone, the stress combination must represent a yield 

condition. Outside the plastic zone, the stress combination must be 

less than that required for yielding. 

. Within the plastic zone, the displacement rate pattern must be related 
to the stresses by the appropriate plastic stress-strain-rate rule. Dis- 
placements at all boundaries must satisfy the appropriate boundary 
conditions. Outside the plastic zone the displacement pattern must be 
that of a rigid body or obey the laws of elasticity dependent on whether 
a plastic-rigid or a plastic-elastic analysis is contemplated. 

. Along the boundary between the plastic and non-plastic zones, appropri- 
ate continuity of stresses and displacement rates must prevail except 
for permissible discontinuities later discussed. 


b. 


Development of the precise solution as described above is considered to be 
too involved mathematically to permit complete solution. By the methods of 
Limit Analysis for a rigid-plastic material, we obtain in effect partial so- 
lutions representing lower and upper bounds for the collapse load. Thus so- 
lutions satisfying subparagraphs a, b and d above define the lower bounds and 
solutions satisfying subparagraphs c and d above define the upper bounds. 
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Yield Condition 


Stress Yield Condition.—The yield condition defines the criteria under 
which plastic flow is possible. The condition can be expressed as an equality- 
inequality between a function of the stress components on the one hand and a 
material constant on the other. Only when the equality exists between this 
function and the constant is plastic flow feasible; when the inequality applies, 
no plastic flow is permitted. In this paper a linear yield condition which ap- 
proximates both the Mises and Tresca conditions is used. This linear con- 
dition is developed from an approximation to the Mises condition for plane 
stress shown in Fig. 2. 

Yield Condition in Terms of Stress Resultants.—As discussed below under 
Equilibrium, the four controlling stress resultants which are considered here- 
in are Nx, Ng, S, and M or in dimensionless form ny, ng, s, and m. From 
the stress yield condition noted in Fig. 2, the corresponding yield condition in 
terms of stress resultants is obtained assuming the stresses in each half 
thickness of the shell to be constant across the half thickness as was done by 
Hodge; (6) the stress resultant yield condition is 


These yield conditions involve four variables. However, the shear stress 
resultant, s, does not combine in any of the inequalities with the other re- 
sultants, and it can conveniently be considered separately. Thus, assuming 
that the shear inequality is satisfied, the remaining yield conditions involve 
three variables. The equalities (b-f) can be pictured in a 3-dimensional space 
as defining a yield surface with all non-yielding combinations of stress re- 
sultants falling inside the surface and points on the surface representing yield 
conditions. The yield surface consists of ten bounding plane surfaces, and 
these are pictured in Fig. 3. 

Because of the uniform loading and shape, the antisymmetric stress re- 
sultants must be zero at the corresponding planes of symmetry. Hence 


at 920, #0 (f), (eg) , and 


at y #0, 8 =0 (h) (22) 


At any transverse section, y = y, the total external longitudinal bending 
moment on the shell must equal the total internal longitudinal bending moment, 


My, resulting from summing N, over the complete transverse section. This 
condition in non-dimensional form is 


aty*y, My 2 2q3KTPL? - y?) sin % (1). (22) 


f 
-lenc +l (c), » (4), 
-len-m-ny c+) (e), -Lomem- ny < +1 (f) (17) 
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If we consider the five stress resultants involved in the equilibrium ofa 
small element as indicated in Fig. 1, it is evident that certain discontinuities 
are permissible while enforcing equilibrium requirements. Thus, in the y- 


direction, discontinuities are permissible for m, Og” and Ng, but ny and s 
must be continuous. For variation in the ¢-direction, discontinuity of n, is 


permissible but m, oo, ng, and s must be continuous. 


Displacement and Strain Rates 


The  pgeseces for rate of internal work per unit area of the shell middle 
surface\9) is expressed in terms of dimensionless stress resultants by 


HX ] (24) 

The strain rate resultants corresponding to the dimensionless stress re- 
sultants are expressed in terms of dimensionless variables by 


Da ™ +n 


au 
= s « 
- wv _ 1 av 


In cases of plastic flow, the plastic action may be confined to a very 
narrow belt. As this belt shrinks to a curve, discontinuities in certain rates 
of displacement and strain may result. Because a discontinuity in W would 
correspond to infinite radial shearing strain in contradiction to the previous 
assumption that the radial shear does negligible work, such a discontinuity is 
W is not allowable. However, discontinuities in U and V and in the first de- 
rivatives of U, V, and W are permissible, provided that due account is taken 
of the energy dissipated at such discontinuities. This may be done by regard- 
ing the discontinuity as the limiting case of a continuous field, as will be illus- 
trated in Section IV. 

The physical conditions inherent in the simple support of the ends of the 
shell impose the requirement that 


at y = +1, V=20, (27) 


Flow Rule 


We consider the plastic strain rate vector defined by 


The Plastic Potential Flow Law then requires that during the plastic flow 
associated with a yield state of stress, the strain rate vector must be in the 
outward normal direction to the yield surface at a point on the surface corre- 
sponding to the given yield state of stress.(11,12,13,14,15) As was indicated 
previously, the yield surface is defined herein in terms of stress resultants 
by Eq. (17) and by Fig. 3. With the yield surfaces consisting of intersecting 
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planes, the flow rule states that if the stress point is on any one plane, then 
the strain rate vector must be outwardly normal to that plane. If the stress 
point is located at the intersection of two or more planes, the direction of the 
strain rate vector must lie between the outward normals to the several inter- 
secting planes. The directions of the strain rate vectors associated with the 
various surfaces of the yield polyhedron are given in Table 1. 


Limit Theorems 


A complete development of the theoretical basis for Limit Analysis is given 
in Reference 9. A synopsis of the discussion applicable to this paper is given 
herein for purposes of continuity. 

A statically admissible stress field is defined for the stress resultants of 
the shell as one which satisfies everywhere in the shell (a) the equations of 
equilibrium (21), (b) the stress resultant boundary conditions (22), and (c) the 
yield conditions (17). As previously noted certain discontinuities in the stress 
resultants are permissible. 

A kinematically admissible velocity field, v*, is defined for the dis- 
placement and strain rates of the shell as one which satisfies everywhere in 
the shell (a) the condition of incompressibility (b) the physical limitations on 
displacements involving continuity requirements and end constraints and (c) 
the condition that positive work must result from the external force, due to 
the radial pressure, acting through the displacement field when considering 
the shell as a whole. As previously noted, certain discontinuities in the dis- 
placement field are permissible. Also there corresponds to each kinematical- 
ly admissible velocity field a plastic stress resultant field, n*, associated 
with the former field by the flow rule. It is noted that although n* is not neces- 
sarily unique, the work done by n* in moving through the strain resultants of 
v* is unique. Further, for each assumed kinematically admissible velocity 
there can be determined a value of the external radial pressure, p**, for which 
the work done by these external forces in moving through the assumed dis- 
placements is equal to the internal work by the stress resultant field n* in 
moving through the strain resultants of v*. 

Lower Bound Criterion.—Any value of the radial pressure on the shell for 
which a statically admissible stress field exists represents a lower bound for 
the collapse pressure, p*. 

Upper Bound Criterion.—Any value of the radial pressure, p**, associated 
as previously described with a kinematically admissible velocity field repre- 
sents an upper bound for the collapse pressure. 


Statically Admissible Stress Fields and Lower Bounds 


Statically Admissible Stress Field I 


It is known that plastic action in ordinary beams occurs with the develop- 
ment of two zones of constant longitudinal stress, one on either side of the 
neutral axis. Since long shells can be expected to act principally as beams, 
in Stress Field I, the half shell 0 < ? < %o is considered in two longitudinal 


zones with the longitudinal stress resultant n, independent of ¢ in each zone 
(Fig. 4). Thus we take 
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TABLE 1: YIELD POLYHEDRON 

Equation Strain Rate Vector 
(€,, , 2 

My -Dg- mel A (1, -1, 0) 


“ng - m= 1 (0, =1, -1, 0) 


“n,+Ng- m= x» (-1, 1, -1, 0) 


ny - = A(1, -1, 

1, 1, 0) 
A (0, 1, -1, 0) 
(0, 1, 1, 0) 
A (1, 0, 0, 0) 
A (-1, 0, 0, 0) 
A (0, 0, 0, 1) 


(0, 0, 0, =1) 
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PLAN OF SHELL STRESS DISTRIBUTION 
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zone A, B ny 


With these assumptions for nx, the equilibrium equations are then solved 
subject to the stress boundary conditions and to continuity requirements on 
the stress resultants along @ = 8. The solutions for the stress resultants are 
for zone A, 0<¢ <8, 


Pl%o B)(sin = y?) 


K 


x 


2 


at 
m = - — 1% sin(% - B)cos 9* 


(sin 9% 9 7-9, ( sinF-sinB) 


Dy = - 1% sin - B) cos - (sin - sin 8)) , (c) 
 3P 


and for zone B, 8 <9<9,, 


Pa(sin %)(1 - x?) 


(a) 
ar? K 


P 8 sin 
ma-77 1% Sin 811 cos (9 %)J - == 


Dg = - x [1+ cos (9 , (c) 


13 PB(sin 9)(@ - 
(a) 


where 


sin 8 - 6 sin 
It will be observed that m and ng depend only upon ¢, whereas nx is a function ~ A 
of y only; the shear stress s depends upon both @ and y. 4 _ 
To permit application of the yield condition, information is required as to 4 
the location of the extrema of the stress resultants, individually and in combi- : 
nation. A summary of the pertinent characteristics for Stress Field I is given ~ * 


below for the typical quarter shell 0 < < 0<y 

(1) ny is invariant with ¢ in each zone. In the y-direction, n, has its ex- 
trema at y = 0. It is negative in zone A and positive in zone B. 

(2) s is negative or zero throughout. In the ¢-direction, its absolute value 
increases linearly through zones A and B from zero to a maximum at 
¢@ = 8. Inthe y-direction, s had its extremum at y = 1. 

(3) m is negative or zero throughout and is independent of y. The absolute 
value of m is a maximum at ¢ = 0 and this absolute value decreases 
monotonically to the edge ¢= ¢. This applies for 0< ¢p <7. 

(4) ng is negative or zero throughout and is independent of y. Its absolute 
value is a maximum at ¢ = 0 and this absolute value decreases mono- 

tonically to the edge ¢ = ¢,. This applies for 0 < $ < 7. 
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(5) It turns out that the combination m - ng needs to be considered in Zone 
A only. The largest value of m - ngin the negative sense occurs at 
¢@ = 0 provided that 8, go, and h comply with 


% Sin B 
B sing —(1-h) (42a) 


In the following we will consider only those values of 8 which satisfy Inequali- 
ty (42a). The above applies for 0< ¢, < %. 


The general yield conditions (17) can now be examined on the basis of the 
characteristics of Stress Fields I. With these characteristics, analysis of 


Fig. 3 shows that all twelve of the yield conditions are satisfied provided the 
following inequalities hold: 


Sly, 9) =-s<1 


(60) 


for n,>0, =ny-m- <1 


for n, <0, G( y,9) + 


Hy,9) >-m-n <1 


(62) 


Let us denote the maximum values of each of these functions by the subscript 
C or A depending upon whether the maximum value occurs at @ = B or at 

¢ = 0. In view of the preceding remarks (1)-(5), the locations of these maxi- 
ma are easily determined. Therefore, the yield condition will be satisfied 
provided that Inequality (42a) and all of 


So = S(l, 8) - 3(1,8) 


Fo = F(0, 8) = - m(8) - ng , 


= G(0, 0) =- Dx,(0) - - 


= H(O) 2 - m(O) - Ng (0) 


(63) 
= N(O) = = ny (0) 


are satisfied. Therefore, in order to determine the best lower bound, we de- 
termine 8 so as to maximize P satisfying (42a) and (63). 

The procedure from here on depends upon the particular values of the pa- 
rameters. If numerical values are assigned to r, h, and $>, then 8 and P can 
be determined so that at least two of Inequalities (42a) and (63) are equalities, 
and the rest are all satisfied. It turns out that Fc = 1 will usually be one of 
the equalities, and that the other one is either Sc = 1, Ga = 1, or Ha = 1. 
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Once a lower bound has been established, the next step is to determine an 
upper bound, using the methods of Section IV. If these bounds are reasonably 
close together, the problem is regarded as satisfactorily solved. However, if 
the bounds are far apart, improvements must be made in one or both of them. 
Generally Stress Field I yielded satisfactory lower bounds for long shells but 
these bounds become unsatisfactory as shorter shells were considered. In 


such cases, it was found that Stress Field I could be generalized as indicated 
in the sub-sections which follow. 


Statically Admissible Stress Field I 


For certain short, deep beams, the shearing stresses tend to become criti- 
cal under increased loading. In Stress Field II shown in Fig. 4, the half shell 
O<¢o< ¢%, is considered in three longitudinal zones with the shear resultant 


in the middle zone taken independent of ¢ and the shear resultants in the other 
zones linear ing. Thus we take 


zone A, qly) , (a) (32) 


zone B, aly) (bd) 


y 
zone C, Ao, qly). (c) 


With these assumptions for s, the equilibrium equations are then solved 
subject to the stress boundary conditions and to continuity requirements on 


the stress resultants along ¢ = a and @ = 8. The solutions for the stress re- 
sultants are for zone A, O< @< a, 


RO - y?) 

R_), R_cos(*-9) _ cosa 


R R cos(P-% Ro 
no PL (2 + cos(u -9) a 1), 
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_ 
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and for zone C, B <%< Po» 


-8) 

RIG - 9)° (e) 

[eos (9, - 9) - 1) * 

h Po ° ( 36) 
ny ° [cos - 9%) , (a) 
(a) 


= 


a( B) sin Po 
(% - 8) sina - a(sin Yo sin 8) 
If a is set equal to @, Stress Field II reduces to Stress Field I. 
Similarly to the discussion of Stress Field I, the pertinent characteristics 


of Stress Field II must now be considered. A summary of the characteristics 
of Stress Field II follows. 


(1) n, is invariant with ¢ in each zone. In the y-direction, nx has its ex- 
trema at y = 0. It is negative in zone A and positive in zone C. 

(2) s is negative or zero throughout. In the ¢ -direction, its absolute value 
increases linearly through zones A and C from zero to the constant extremum 
over the entire interior zone B, a<¢<f. In the y-direction, s has its ex- 
tremum at y = 1. 

(3) m is negative or zero throughout and is independent of y. The absolute 
value of m is a maximum at ¢ = 0 and this absolute value decreases montoni- 
cally to the edge @ = ¢,. This applies for o< $9 < 7. 

(4) ng is negative or zero throughout and is independent of y. Its absolute 
value is a maximum at ¢ = 0 and this absolute value decreases monotonically 
to the edge ¢ = ,. This applies for o< ¢o < 7. 

(5) The combination m - ng needs to be considered in Zone A only. The 


largest value of m - ng in the negative sense occurs at @ = 0 provided that 
a, and h comply with 


~ ein (9, - 2) 105 


With these stress resultant characteristics, the arguments previously ad- 
vanced in the simplification of the yield condition for Stress Field I apply 
equally to Stress Field II except that Inequality (42) must be enforced in lieu 
of Inequality (42a). Hence Inequalities (63) replace the original yield condition, 
and we seek to determine both a and f so as to maximize P while satisfying 
(42) and (63). The best lower bound on P is thus obtained. 

As previously the procedure now depends on the particular values of the 
parameters. After numerical values are assigned to r, h, and g,, aand B 
and consequently P can theoretically be evaluated so that at least three of the 
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Inequalities (42) and (63) are equalities and the rest are all satisfied. How- 
ever, because of the rather complicated equations involved this method did not 
prove practical. Instead lower bounds were obtained by trial, assigning values 
to a and B on the basis of previous results and using the best resultant lower 
bound, 

Again as previously, the lower bounds obtained by this method were judged 
to be satisfactory only if close enough to an upper bound determined by the 
methods of Section IV. In the procedure adopted for computing lower bounds, 
Stress Field II was tried only in those cases where the results of a previous 
trial using Stress Field I indicated Sc to be the controlling stress resultant. 
This was true of the shortest and thickest shells considered. In these cases 
the greater shear capacity of Stress Field II resulted in a higher and satis- 
factory lower bound. [If the lower bound from Field I was judged to be too low 
and the controlling stress resultant combination was not Sc, Stress Field III, 
representing a further generalization of Stress Fields I and II was tried. This 
stress field is taken up in the succeeding sub-section. 


Statically Admissible Stress Fields II 


Stress Fields II] are derivative from Stress Fields I and II previously con- 
sidered which specified the stress resultants for the shell of subtended angle 
2¢%5 with symmetry about the center line @= 0. Now taking a shell of subtend- 
ed angle 2n§, with n an integer we consider the shell as consisting of n ad- 
joining sections, each of subtended angle 26,, which extend from the near 
edge $= n@, to the far edge @ = -n@9. We assume that the stress field in 
each such section of subtended angle 26, is either Stress Field I or II as ob- 
tained for a shell of subtended angle 26,. The near edge ¢ = n@, will super- 
impose on the edge of a narrower shell Ba 6, and similarly the far edge 

= - n@, will superimpose on the edge ofa aarrower shell g-= = -6. Since 
the boundary conditions at ¢ = *n@, are the same as at = - 74 for the 
narrower shell, the boundary conditions at d= *n@, for the full. shell are 
satisfied by the assumed stress fields. At all internal junctions between 
sections, ¢ = (n - 2K)@o, with K also an integer, the stress resultants, s, ng, 


om 
m, and od are each zero and thus necessary continuity conditions are com- 


plied with. Hence, subject to compliance with the yield condition, Stress 
Fields I and II developed for a shell of subtended angle 26,, are statically ad- 
missible for shells of subtended angles 2n@,, n integral. Thus with the use of 
Inequalities (42), (42a) and (63), lower bounds can be developed as described 
in the preceding sub-sections. As previously indicated Stress Field II was 
considered only after previous use of Stress Field I gave an unsatisfactory 
lower bound with a controlling stress resultant combination other than Sc. In 
all such cases Stress Field III yielded a satisfactory lower bound. In general 


this stress field was found to give the controlling lower bound for the shorter, 
thinner shells. 


Kinematically Admissible Velocity Fields and Upper Bounds 


Kinematically Admissible Velocity Field I 


The displacement pattern of this field is similar to that associated with 
collapse of a simple beam. As shown in Fig. 5, a plastic belt between y = “yj 
and y = yy is assumed with the remainder of the shell rigid. Displacement of 
each rigid section is taken to correspond to rotation about a horizontal axis in 
each end passing through the shell arc at ¢ = : 8. With displacements 


7 
i 
. 
5 
ay 
x 
= 
of 
od 
¥ 
if 
E 
we’ 
a 
ag 
* 


1706-16 EM 3 July, 1958 


NOIL93S 
1% =f 40 
Saxy 


NOIL93S 


KINEMATICALLY ADMISSIBLE VELOCITY FIELOT 


PL-=Plastic Line 
P.C=Plastic Curve 


NOIL9D3S 


KINEMATICALLY ADMISSIBLE VELOCITY FIELO TI 


FIGURE 5 


| 
8 | 
> PLAN > 
KINEMATICALLY ADMISSIBLE VELOCITY FIELD I se 

° é 

| R, Re a| R, 

= PL | Po PL 
\e 
b=. 
\ 
| 


CIRCULAR SHELL ROOFS 1706-17 


symmetrical about ¢= 0 and y = 0 the quarter shell 0< ¢ < ¢, O< y< lis 
typical and the expressions for dimensionless displacements are 
in R, < y 1 


U ar(cos B- 9), 


a 


a 
and in P, O<y “ = cos 


ar(cosR- cos g) & 


a Y1 
= (sin 9) (1 - > - 1 ) 


(44) 
It is noted that these displacements obey the physical constraints of the simple 
supports with V = W = 0 at y = 1 and that W is continuous throughout the shell 
as required. The accompanying strains are derived from the displacements 
and the rigid body nature of the displacements in R is verified. Thus, in R, 


Hx = 0 . 
However, in P, 


= (cos cos 9), 


for 


for % . 


With the above combinations of strains, the applicable plastic regimes for 
P are for 0< ny = -1 and for BS. go, ny = 1. These regimes corre- 
spond to faces K and J respectively of the yield polyhedron, Fig. 3. 

The internal work done in the quarter shell is designated I' and we define 
I by 


2/3 TK 


I (n, Dg + mHX)d8, (45) 


We use the previous values of nx and €x in P to obtain 


yz { pr - (cos cos g) Lady age SP (cos B- cos g)L ady ag} 
Yl 


Lear [9 cosB- sin 9, - 2B cos B+ 2 sin 
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The work done by the external forces on the quarter shell is designated E' 


mds . (46) 


For Field I, 


E® Pla { = (cos g)(1- 4 - “Sy, (1-y)avay } 
‘ Plaa(3 - Sin % 
r 
For P to be an upper bound, E = I, from which we obtain 
cos B= sin - 2B cos + 2 sin 8) 
(3 - yy") sin % 
Minimum values of P with respect to y, and B result for yy = 0 and & = 


do 


These values of y; and # are used to give the expression for an upper bound 
on P, 


9 
2r?(2 sin - sin %) 
sin ?, 


(47) 
It is noted that by taking y, = 0, the plastic belt has been shrunk to a curve, 
This results in a discontinuity in U which is permissible. Although €, be- 
comes infinite, the internal work remains finite. 

If the upper bounds obtained from (47) are reasonably close to the lower 
bounds obtained in Section III for a particular value of the shape parameters, 
no further effort towards a better upper bound is made. In general, Velocity 
Field I gives satisfactory upper bounds for the longer shells considered. How- 
ever for shorter, thinner shells of larger central angle, the departure of the 
upper bound from known satisfactory lower bounds tends to become unreason- 
ably large and additional velocity fields must be investigated. Such fields are 
taken up in the succeeding subsections. 


Kinematically Admissible Velocity Field I 


In the displacement pattern of Field II, the complete shell is divided direct- 
ly into 6 rigid sections separated by plastic curves as shown in Fig. 5. Region 
Rg rotates as a rigid body about a plastic line at @ = 0 and undergoes a simul- 
taneous tangential displacement. Region Rj rotates as a rigid body about an 
axis at the end y = 1 which is parallel to the normal to the shell at ¢ = 0 but 
which is displaced laterally therefrom. In addition region R, undergoes a 
rigid body twisting motion. Although the chosen displacement field was set 
up originally on the above physical basis it suffices now to verify the rigid 
body nature of the displacement pattern and the satisfaction of continuity and 
boundary requirements. Again we consider the quarter shell OK oS %o and 
0<yS1. The displacement pattern in Ro, 0X  S do, OS y yq is chosen 
as 


TK 
Then 


CIRCULAR SHELL ROOFS 


VEZ (a- cong), sin ?, (48) 


and it is evident that € , = Ed = Y = HX = 0 in Rg as required for a rigid body 
motion. The displacement pattern in Ry, < y S1 is chosen as 


(Ap - sin + B) 
1 


= - y) 
(A cos 9) 


a(l - y) 

The internal strains with this displacement field are all zero. With regard to 
W, it is noted that W is continuous in each of the regions R, and Rpg and at all 
junctions between regions as required. At y = 1, the physical constraints are 
satisfied as V = W = 0. 

With the displacement pattern of Field II, internal work is done along the 
plastic curves. The discontinuities along these curves are investigated for 
consequent internal work as noted below. 


Discontinuities Along $= 0, 0< y < yy.—Pertinent discontinuities applic- 
able to the quarter shell are* 


I,, the internal work expression (45) corresponding to these discontinuities, 
is minimized by setting 


(50) 
Then, 
La a hy) 
1“ ° 


Discontinuities Along ¢ = 0, y1< y < 1.—Pertinent discontinuities applic- 
able to the quarter shell are 


a(l - A-1 
Vv 
Ig, the internal work expression corresponding to these discontinuities is 
again minimized by evaluating A as in (50) and we obtain 
laah(1 = 
2 ~ +h) 
Discontinuity Along y = yy, 0< @< $.—A discontinuity in U exists, name- 
ly 
- sing + B). 


The value of B is selected so as to minimize Iq, the internal work expression 
corresponding to the discontinuity in U. The form of the expression I is 
dependent on the value of ¢,. Thus, for 


*The symbol denotes a discontinuity in the quantity immediately preceding. 
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2 
- laar %o Po 


2 
laar Fo Fo 
{1 - 2 cos +008 %* J 
The expressions defining 7, 4, $1, are 


Y 
+h 


H< 26, = - 


-siny=0O , 


(51) 


Fo 


%- 


- sin sin 20, 


54 
26+ wind 0. (54) 


The expression for the external work associated with Velocity Field II is 


Plas 


1 
singag ty, sing ag}. 
1 
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For P to be an upper bound, E =I, + Ig +13. We can write P in terms of 
y, and three constants which are independent of y;. Thus 


Cc 
. 
- 


Cc 
1 
Then, P will have its lowest value with respect to y; when y; = 0. The re- 
sultant expressions for the upper bound P** are 


l- cos 9 
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As noted in the preceding subsection, Velocity Field II has been investigat- 
ed in an effort to obtain better upper bounds than those yielded by Field I par- 
ticularly in the area of the shorter, thinner shells of larger central angle. It 
is found that Field II does give better upper bounds than Field I in this region 
of the shape parameters and in an isolated instance does give the controlling 
value of the upper bound. However, generally the upper bounds from Field II 
are still not considered satisfactorily close to known lower bounds, and it is 


necessary to investigate Velocity Field III, a generalization of Field II taken 
up in the subsection directly following. 


if 
= 
where A = 
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Kinematically Admissible Velocity Field III 


This field is derivative from Velocity Field II in which it was found that 
minimizing of the resultant upper bound enforced the shrinkage of a central 
rigid section -y)< y < y, to a curve at y = 0. In Field II, the shell is divided 
into 5 rigid sections separated by plastic curves as shown in Fig. 5. However, 
the rigid section Rg, -u << uw, is completely motionless and only the 
sections Rj, andu< S po, are in motion. Plastic lines occur 
at @ = *w and plastic curves at y = Oin the ranges -Q@<@ -u and 
USS $%. We consider the quarter shell in the first quadrant and select for 
the displacement pattern in Ro, w and y< 1, 


(56) 
and in Ry, u andO<y<1, 


U sar [Al@- yp) - sin * BJ 
y LA - cos (? 


a(l y) 


sin (9 - 


Investigation readily reveals that the derived strains are zero in each of 
regions Ry and Rg, that V and W vanish at the ends, and that W is continuous 
throughout. 

We now set ¢ - u equal to 6 and gp - HM equal to 6, and compare Fields I 
and III using y; = 0 in Field II. It is apparent that the expressions for the dis- 
continuities and the external and internal work in Field II are identical with 
those previously obtained for Field II except 6 and 6, replace ¢ and gp. Ex- 


pressions can therefore be written for the upper bound P** analogous to (55) 
as 
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By the use of equations (58), P can be minimized with respect to 6, for par- 
ticular values of the parameters r and h and the minimizing values of 6, es- 
tablished. By this procedure, collapse modes involving collapse of only 

portions of the shell are developed. 


Solutions for Bounds 


Computations were carried out to give upper and lower bounds for the col- 
lapse pressure parameter P with a wide range of values for the parameters 
%o, r, andh. With ¢, taken first as 30° and then as 90° and h taken in turn 
equal to .02, .005, and .00125, r was varied from 0.2 to 4.0 in 5 steps. The 
results are plotted as 6 pairs of curves in Figs. 6 to 11. Tabulation of the re- 
sults with associated pertinent developmental data is given in Tables 2-7. 
With regard to Velocity Field II, the defining ranges of 6,, 2” < , 
26 < 6,5 27, 0K 69 < 26, have been designated Cases 1, 2, and 3 for pur- 


poses of tabulation. Special data developed for Field III are tabulated in 
Table 8. 


CONCLUSIONS 


1. For the case of a simply supported cylindrical shell without edge beams, 
collapse pressures and collapse modes have been obtained. 


a. By the use of methods of Limit Analysis for a plastic-rigid structure, 
lower and upper bounds for the collapse pressure parameter, P*, 
and P**, have been developed. Subject to the validity of the basic 


assumptions, if P = ; (P* + P**) is adopted as the solution for the 


collapse pressure parameter, the results indicate this solution to be 
accurate to within - 2 to 25 per cent for the wide range of the shape 
parameters, $,, r, h considered. 

b. It is realized that the results obtained herein for P* and P** do not 
define precisely the collapse modes. However, review of the con- 
trolling stress and velocity fields as tabulated in Tables 2-7 indi- 
cates that when the controlling lower bound is defined by Stress 
Field I the controlling upper bound is usually defined by Velocity 
Field I and that when the controlling lower bound is defined by Stress 
Field III, the controlling upper bound is usually defined by Velocity 
Field III. Now Stress Field I and Velocity Field I have been based on 
characteristics associated with the collapse of simple beams. Since 
the results for lower and upper bounds obtained from these fields 
for long shells are in close proximity it seems reasonable to con- 
clude that collapse of long shells will approximate that of simple 
beams. Similarly Stress and Velocity Fields III each appear to be 
associated with partial collapse of the shell along the edges. Hence 
Velocity Field III is considered to be a reasonably close approxi- 
mation to the collapse mode of those shells for which the controlling 
value of the upper bound P** was obtained from Field III. 


2. Considered assumptions made in the development of this paper are now 


reconsidered to assess their possible effect on an extension of the results to 
real shells. 
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TABLE 2: % = 30°, h = .02 


LOWER BOUNDS 
Tr ol 
011 
S.A. Meld I I I 
Bue Criterion Fo s Gy, Fo Ga Fo Gy 
B 15° 15° 14° 30! 
a - 


UPPER BOUNDS 


r 1.0 2.0 
p* 200288 20115 2072 2236 


K. A. Field z I I I III 
Case - - - 2 


- - - - 30° 
16° 


3° 
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TABLE 3: go = 30°, h = .005 


LOWER BOUNDS 
8 1.0 2.0 4.0 
P 200275 20103 20339 20477 -097 +153 
S. A. Field I I I b I III 
Bea Criterion Fo= Gp Fo Gp Fo # Gp Fc = Gp Hy So 
B 15° 15° 16° 17° 240401 7°30" 
a - - - 
UPPER BOUNDS 
r 2 4 8 1.0 2.0 4.0 
-00288  =.0115 072 .134 +235 
Ke Ae Field I I I I III III 
Case - - 1 2 
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TABLE 4: 30° h = .00125 


LOWER BOUNDS 


1.0 2.0 
-0210 +0250 2047 

I I III 
Fo Ga Fo Fo Gy 
23°20" 25°20' 8°29" 


15°00" 
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.4 1.0 2.0 
15 +0300 -0730 
I III III III 
- 1 1 1 
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TABLE 5: 9, = 90° h #® .02 


LOWER BOUNDS 
p* 0980 193 2272 
S. A. Field I I III III III III 
Criterion Fo #G, Fos Ga Fo SG, Fo Fox Ha 
8 54°40! 76°00! 22°30" 23°10" 20° 27° 
9 45° 45° 300 30° 
° 
UPPER BOUNDS 
4 1.0 2.0 4.0 
K. A. Field I III III III III III 
Case - 1 1 1 2 2 
% - 64°10" 45020" 41°10" 30° 24° 
$ - - - - 30 4°50" 
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hb = .005 


TABLE 6: 9, = 90 
LOWER BOUNDS 
4 8 1.0 2.0 4.0 
III Ill Ill III 
26°10! 35030! 39°10 24°40! 


UPPER BOUNDS 
4.0 


r 24 1.0 
p** 0164 0310 0591 0724 .134 
Ke Ae Field III III Ill III III III 
Case a 1 1 1 a 2 

63° 45°40! 33°30" 28°40" 20°30" 16°00" 
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TABLE 7: 9 = 90° h # .00125 


LOWER BOUNDS 
Tr e2 1.0 2.0 4.0 
P* 200469 =.0103 20250 +047 095 
S. A. Field III III III III Il Ill 
Bea Criterion FotG, Fo Fot Gg Fo Hy 
26°56" 37°04" 23920! 250201 8°29" 12°428 
9 45° 5° 30° 30° 15° 15° 
UPPER BOUNDS 
r 1.0 2.0 4.0 
.00791 .0300 0372 -0730 
K. A. Field III III III III III III 
Case 1 1 1 1 1 1 
% 45°30" 33°40" 23°00" 21900" «13920 100308 
2 - - - - - 
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TABLE 8 
Date for Velocity Field III 
b A Y 6 Case 1 Case 2 Case 3 
202 980392 19°50" 6° 12% 
005 995025 9°50% 3°45! 19°40*<G, 730" 


00125 «998752 4958" 1957.95 99569<8, 
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Chairman of the dissertation Guidance Committee for his many helpful 
suggestions during the development of the dissertation. 
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As we have dealt with a rigid-perfectly-plastic material, we have 
assumed no deflections prior to the onset of plastic action. While 
this approximation may be close for short shells the elastic de- 
flection of long shells up to the onset of plastic flow may be so large 
as to make the shell non-usable for its intended purpose. The as- 
sumption of perfect plasticity which assumes no strain hardening is 
not believed to offer serious objections for ductile materials such 
as soft steel which undergo substantial plastic flow prior to onset of 
strain hardening. 

The yield condition used herein represents a linear approximation to 
the Mises yield condition. However, since the Mises condition also 
represents the behavior of actual materials only approximately, it 
does not appear unreasonable to use this linear yield condition for 
ductile metals. The linear condition used is closest to the Mises 
condition where the shear is small and its validity is therefore be- 
lieved to be greatest for small shear. Of course the linear condition 
can be readily revised to approximate the Mises condition more 
closely for larger shears, but this would complicate the work of at- 
taining solutions. 

In the development of equilibrium equations, certain stress re- 
sultants were neglected. For very short shells, the omission of the 
longitudinal bending moment from equilibrium considerations ap- 
pears to require critical evaluation. Thus very short shells of small 
curvature could be expected to act as plates supported on the ends. 
However the effects of shortness and curvature tend to magnify and 
lessen respectively the importance of this longitudinal bending 


moment and no conclusions are drawn at this time as to the effect of 
this omission. 


. The expressions used for strains and the equilibrium equations pre- 


supposed small deflections. Although this is not in accord with the 
deflections associated with collapse no conclusions are drawn as to 
the resultant effect on the accuracy of the solutions. 


. With regard to the applicability of the results obtained herein to 


Shells built of reinforced concrete, it is considered that present 
information on yielding of this material under the three dimensional 
stress system of the shell is insufficient to permit assessment of 

the extent of the approximation involved. However, it is believed 
that with the development of necessary yield criteria for reinforced 
concrete, the methods of limit analysis as used herein can be utilized 
to advantage to develop the desired data on plastic collapse. 
Extension of the methods of this article to cover other end and edge 
conditions of the shell is suggested. 
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List of Symbols 


radius of shell 


t 
dimensionless thickness-radius parameter = 2a 


circumferential bending moment in non-dimensional form 


any integer 


longitudinal direct stress resultant in non-dimensional form 


= circumferential direct stress resultant in non-dimensional form 


uniform radial pressure on shell 


lower and upper bounds on p 


dimensionless radius-length parameter 


membrane shearing stress resultant in non-dimensional form 


thickness of shell 


displacements in longitudinal, circumferential, and radial directions 


longitudinal coordinate 
longitudinal coordinate in non-dimensional form 


radial coordinate 

constants defining Velocity Fields 

total internal plastic work in shell 

reduced expression for work done by external forces 


displacement vector 


Ha, Na, Sc = critical combinations of stress resultants relative to 
yield condition 


K = functions of ,,@, 8 defining the stress resultants 


distance between centers of faces of the equivalent ideal sandwich 
shell =t 


reduced expression for internal plastic work 
constant of Mises yield condition 

length of shell 

circumferential bending moment per unit length 


limit bending moment per unit length = 2 V3THK = V3 t2K 


= membrane direct stress resultants per unit length 


we 
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No = limit direct stress resultant per unit length = 23 TK = 273 tK 


P= TY3RT = dimensionless pressure parameter 


P = value of P at collapse 
P*, P** = lower and upper bounds of P 
Q = radial shearing force per unit length 
= membrane shearing force per unit length 
So = limit membrane shearing force per unit length = 2TK = 2tK 
zy = thickness of one face of the equivalent ideal sandwich shell = t P 
U == = longitudinal displacement in non-dimensional form 
Va : = circumferential displacement in non-dimensional form 
We = = radial displacement in non-dimensional form. 
a = 


= constant associated with rigid body rotations in Velocity Fields 
a, B,H = values of ¢ defining boundaries between longitudinal zones 

Y, 6, 4, $> = values of g, used with Velocity Field I 

Y, 6, 01, 92 = values of ¢, used with Velocity Field I 


> = angular coordinate in circumferential direction 
25 = central angle subtended by the shell arc 

6 =@-u 

Ex, €g = unit elongations 


Y= "xo = shearing strain 

X=Xqg = rate of change of curvature in circumferential direction 
Ox, %, G2 = normal components of stress 

Tx Tz Tzx = shearing components of stress 


=Tx¢ 
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EFFECTIVE WIDTH OF THIN RECTANGULAR PLATES 


Morris Ojalvo,* A.M. ASCE and Frederick H. Hull** 
(Proc. Paper 1718) 


SUMMARY 


The behavior of simply supported rectangular plates subjected to edge 
compression in one direction is investigated by means of load-shortening 
tests. The data are presented by non-dimensional curves and compared with 
the theoretic results of Levy(1) and Coan(2) and the experimental data of 
Sechler,(3) Winter,(4) and Lahde and Wagner. (5) 


INTRODUCTION 
Edge supported rectangular plates are subject to buckling when com- 


pressed by edge forces in the middle plane of the plate. The theoretic buckling 
stress when the load is parallel to a pair of opposite edges is given by: 


6,F 


2 
Yen = 


The constant K depends upon the length ratio of the sides of the plate and the 
degree of rotational restraint at the edges. In actual tests buckles develop at 
stresses smaller than those given by equation (1) because of the deviations 
from flatness found in all real plates. 

The buckling of a plate does not signify the limit of its load carrying 
capacity. After buckling the plate continues to carry increased loads but at 
reduced stiffness. That is, for any additional load increment the plate 
shortens at an increased rate. Meanwhile the average stress distribution 
over the thickness of the plate becomes non-uniform on any section 
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transverse to the load with the larger compressive stresses at the edges of 
the plate. This has given rise to the effective width concept.(7) A useful 
definition for the effective width is given by equation (2): 


— 
e tee 
P is the load carried by the plate at a unit shortening of e. 
The effective width is always smaller than the total width, b. Its signifi- 


1 
cance may be seen if equation (2) is multiplied by b and rearranged: 


PAt _ Tove Je _ be (3) 
® 


The ratio be to b equals the ratio of the secant modulus of the average stress- 
unit shortening curve of the plate to the modulus of the plate material. With 
this ratio known, it is possible to evaluate closely the contribution of com- 
pressed plate elements when they are part of semi-monocoque structures or 
occur in members formed of thin sheet metal. Theoretical analysis(1,2) 
shows that the unit longitudinal strain at the edges parallel to the compres- 
sive load is variable. “e” is the average value of the longitudinal unit strain 
over the length of one buckle. It may be evaluated by computing the unit 
shortening for a long plate having several buckles. 

The plates tested were supported along the edges in a manner providing 
little rotational restraint. The longitudinal edges were allowed to move in the 
plane of the plate. 


Description of the Test Procedure 


Each of the two jigs used in the tests consisted of a six-inch I-beam with 
two guide bars (figure 1). Two of the I-beam flange edges, a, were carefully 
machined so that the surfaces were smooth and lay in the same plane. The 
outer surfaces of the I-beam flanges, b, were partly machined to provide a 
smooth surface for the guide bars, A, to bear upon. The surfaces, c, of the 
guide bars were carefully made plane so that the surfaces c and a could form 
the parallel sides of a groove along which the longitudinal edges of the plate 
slide. The bolts, B, clamp the guide bars in position to accommodate the 
different thicknesses of plate tested. A 1/4 inch base plate with a shallow 
slot cut in it keeps the lower edge of the plate in its plane while a 3/8 inch bar 
with a similar slot is used at the top of the plate (figure 2). The clearance 
with the jig at each end of this top bar was made as small as possible, it being 
about 1/32 of an inch for the thinnest (.025 inch) plates. 

The plate material in all of the tests was 24S-T3 aluminum. The mechani- 
cal properties in compression were obtained from stress-strain curves of 
small coupons restrained from buckling. Some of the plate tests were con- 
ducted with the load parallel to the grain and some with load transverse to the 
grain. All of the plates were 5.80 inches wide. They were inserted into the 
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jig as shown in figure 3. Table 1 summarizes the dimensions of the test 
specimens and the mechanical properties as determined from the stress- 
strain curves of the coupon tests. In machining the plates special care was 
taken to make the ends parallel. 

Figure 4 shows one of the long plates being tested. Dial gages reading 
directly to the nearest 1/10000 of an inch are placed in the plane of the plate 
to measure the shortening. A rigid bar, C, is placed between the loading 
head of the testing machine and the slotted top bar to fascilitate this mea- 
surement. In practically all cases the dial readings were in close agreement. 
The average was used to give the shortening. The testing machine was of 
400,000 lb. capacity operating in the 0-16,000 lb. range. The least count on 
the load dial was 20 lbs. 

Consideration was given to the possibility that a substantial portion of the 
load transferred to the jig by means of friction between the longitudinal edges 
of the beam and the grooves of the jig. This would mean that the load sup- 
ported by the plate is greater at the top and that as a consequence the stiff- 
ness of the plate would appear greater than its true value. To minimize this 
effect the longitudinal edges were lubricated with Molykote type G thinned with 
a fine grease. To test the effectiveness of this procedure SR-4 type A-12 
strain gages were attached to the front and back of a .082 inch thick plate at 
the locations shown in figure 5. These were wired to give the longitudinal 
strain in the middle plane of the plate. The proximity of the gages to the 
edges meant that the normal stresses in the transverse direction were small 
and that therefore the strain reading was a good measure of the longitudinal 
normal stresses. Before buckling took place the strain readings indicated 
practically no strain variation along the edges. For higher loads the strains 
in micro-inches are plotted against the load on the plate (figure 5). The maxi- 
mum load of 4800 lbs. applied in this test was about half the load at failure. 
There is no indication that the readings obtained from the bottom gages are 
generally lower. The spread in readings is expected from the theory(2) which 
predicts variations in the longitudinal stress after buckling. The theory pre- 
dicts much higher edge stresses at the center of a buckle as compared with 
edge stresses at a nodal line. The best indication that the frictional transfer 
was not large is found by comparison of the performance of the long (46.40 
inches) plates with the short (23.85 inches) plates. In figures 6, 7, 8, and 9 
where be/b is plotted against e/ecr, the results for the long plates usually 
fall within the band of scatter established for all plates of the same width to 
thickness ratio. 


Observed Behavior of the Loaded Plates 


All of the plates showed buckles at loads smaller than the critical buckling 
load (figure 10). The buckles initially numbered 8 for the long plates and 4 
for the short ones. As the load was further increased the amplitudes of the 
buckles increased and also the lengths of some of the buckles increased while 
others decreased. This behavior was followed by an increase in the number 
of buckles. Sometimes the increase in the number of buckles occurred with 
an audible snap. 

In the case of one .025 inch thick plate, buckles near the edge were noticed 
in addition to the primary ones when the critical strain had been greatly ex- 
ceeded. Figures 11 and 12 are photographs of this kind of buckling as viewed 
from two different angles. 
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Failure was always due to a tearing at the top along longitudinal lines be- 
tween the top bar and the grooves. The thicker plates showed a little of the 
buckled shape remaining throughout the length after removal of the specimen 
from the jig, thus indicating some plastic flow under combined stress at the 
higher values of the load. 


Discussion of the Test Data 


To determine the average stress the load was divided by the area based 
upon the full width (5.80 inches). It was felt that this was preferable to esti- 
mating the portion of the load carried by the edge strips within the grooves 
and basing the average stress on the free width of the plate. The portion of 
the width in the grooves is just under 4% of the entire plate width. 

The critical buckling stress (theoretical) was taken as the buckling stress 
for a square plate simply supported along the edges. In the formula for it 


given below the full width (5.80 inches) was used for b: 


2 


_WeE t (4) 
Cer. = 301-4?) b* 


Poissons ratio was taken as .32. The critical strain is determined by 
dividing the critical stress by the modulus of elasticity. 

The test data are presented in figures 6, 7, 8 and 9 in non-dimensional 
form by plotting be/b against e/ecy,. The solid lines represent an average 
curve while the band of scatter within the dotted lines includes the results of 
5 to 7 load-shortening tests. The results of the tests on the long plates are 
indicated by dot-dash lines within the band of scatter. These are similar to 
the curves obtained from the tests on the shorter plates. 

The wide band of scatter at low values of e/ecr, shows the influence of 
initial deviations of the plate from flatness. Theoretic work considering 
small, regular deflections of square plates(2,8) indicates that this is to be ex- 
pected. Deflection measurements made in the manner shown by figure 13 
were taken on one of the .025 inch plates and on one of the .082 inch plates. 
Maximum deflections were about 60% of the thickness for the .025 inch plate 
and 20% of the thickness for the .082 inch plate. The initial deflections were 
entirely on one side for the .082 inch plate. The .025 inch plate had deflec- 
tions below the plate parallel to one longitudinal edge and above the plate 
parallel to the other edge. 

Figure 14 compares the results of the present investigation (dotted curves) 
with other theoretical and experimental results. The tests of this paper give 
evidence of an additional dependence of the data on the b/t value. This effect 
does not manifest itself in the theoretical solutions and is not indicated in the 
results of other investigations. 

At this point it is profitable to revue the boundary conditions assumed for 
obtaining the theoretical results and the methods used to obtain the experi- 
mental results of figure 14. 

Levy’s solution is for a square plate with the edges free of rotational re- 
straint. All edges are constrained to remain straight, there are no shears 
on the edges, and there is no resultant transverse force on the longitudinal 
boundaries. He has obtained solutions for the displacement normal to the 
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plate in terms of doubly infinite series for a number of different compressive 
loads. These satisfy the boundary conditions and Von Karman’s two differen- 
tial equations for the large deflection of plates:(9) 


axt* oy* ax* ay? 


ax? dy? 


+ (5) 


The results were obtained by assuming all but the most significant six terms 
of the series to be equal to zero. The error due to this is quite insignificant. 

Coan’s solution differs from Levy’s in that the longitudinal edges are not 
constrained to remain straight. This corresponds most closely to the present 
tests. In this solution only three terms of the series are assumed to be dif- 
ferent from zero. This too is believed to introduce only minor errors. 

Sechler’s results were obtained by compressing simply supported plates 
to failure. The effective width was determined by dividing the load at failure 
by the product of the yield strength and the thickness. 

Winter evaluated the effective width by testing beams made of thin gage 
sheet metal. In the elastic range the location of the neutral axis was deter- 
mined from the assumption that the strain is linear over the depth of the beam. 
The effective width of the compression flange was then evaluated from the 
condition that the neutral axis was also the centroidal axis of the section. In 
the plastic range a stress distribution as in a plastic hinge was assumed and 
the effective width was computed from the condition that the resultant of hori- 
zontal forces on the transverse section disappears. The equation that best 
fits the data of both Sechler’s and Winter’s tests is given by: 


(6) 
- -298 


The tests of Lahde and Wagner were made on long plates with the edges 
fixed against rotation. The changes in the plate length in the longitudinal and 
in the transverse direction were measured. The load carried by the plate on 
a transverse section was computed from numerous slope measurements on 
the surface of the plate. Since there are shearing forces on the edges of the 


plate the computations give the average values for the longitudinal load car- 
ried by the plate. 
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RECTANGULAR PLATES 
CONCLUSIONS 


1. Initial deflections of plates having the same width to thickness ratios 
cause large variations in the effective widths at low values of e/e€cr.- As 
e/ecr, becomes larger the variation in the effective width grows smaller. 


2. In general, for any value of e/ecr, the value of be/b will be larger the 
lower the width to thickness ratio of the plate. 


3. The results of Sechler’s and Winter’s tests are conservative for design 
purposes except when b/t is large and e/ecr, is small. 


4. It is desirable that further tests include: 


a) extensive determinations of the initial deflections for the purpose of 


correlating the behavior of the plate with the magnitude and shape of 
the initial deflections; 


b) specimens having the same b/t ratios as in these tests but different 
widths and thicknesses. 


Symbols 


width of the plate normal to the direction of the load. 


effective width of plate. 

t3 
plate constant equal to ina 
modulus of elasticity. 


unit shortening of the plate in the longitudinal direction over a length 
of one buckle. 


Ocr, divided by the modulus of elasticity. 
stress function. 

- constant, 

- compressive load carried by the plate. 

- plate thickness. 


- displacement of the middle surface of a plate perpendicular to its 
plane. 


- coordinate. 

- coordinate. 

- poissons ratio. 

- average compressive stress over entire width of plate, b. 


- the critical stress at buckling (theoretical) for a square flat plate 
simply supported at the edges and under uniaxial edge compression. 
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= cosh Ya + } 
Rm = coth Gb + 
t = TEMPERATURE DIFFERENCE = t(4%)-t(-$) 
{ (340) cosh + | 
W = DEFLECTION 
Way, = + Ran 
Oo = COEFFICIENT OF LINEAR EXPANSION 
Xm = 
= 
= Poisson's RATIO 
+% a coth ta 


ABSTRACT 


The problem of bending of a rectangular plate, elastically supported on 
variable rigidity beams, subjected to a transverse loading or temperature 
gradient is solved in this paper. 

Specific combinations of support are examined and, in addition, a detailed 
solution is given for a square plate with one edge free. 


INTRODUCTION 


The problem of the bending of an elastically supported rectangular plate 
subjected to transverse loading or to temperatures varying linearly through 
the thickness will be considered in this paper. The plate is assumed to be 
supported on beams which in turn are simply supported at their ends. The 
four beams may have variable rigidity along their length, and this rigidity 
may vary from beam to beam. The deflection is determined and from this the 
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moments and stresses can be readily calculated according to the classical 
theory(1) of thin plates. The plane stress problem, associated with the sym- 
metric component of the thermal gradient will not be considered here. 


The method of solution is to consider the superposition of deflections of 
the following: 


1. Rectangular plate, simply supported on all four edges, subjected to the 
given temperature distribution or loading system. 


2. Rectangular plate, subjected to an arbitrary line load on one edge and 
simply supported on the other three edges. 


3. Simply supported beam of variable rigidity, subjected to an arbitrary 
load distribution. 


The satisfaction of equilibrium and compatibility along the plate edges 
yields the unknown line loads, from which the complete solution can be ob- 
tained. 

Solutions have been previously obtained(1, 2, 3, 4, 5) for the plate under 
transverse loading and supported by constant cross section beams arranged 
either in grid like fashion, or paired in such a manner as to have equal 
rigidities on opposite plate edges. This paper represents a generalization of 
the edge supported plate to include thermal gradients and arbitrarily assigned 
variable cross section supports. 

The technique of superposition of various elementary problem types leads 
to an understanding of the nature of the support reaction, which in turn points 
up and allows detailed investigation of certain difficulties inherent in the ther- 
mal problem. In addition this technique is useful in the establishment and 
evaluation of the accuracy of approximations, and in obtaining separate checks 
of the numerical calculations. Alternatively, the boundary conditions may be 
written in terms of the total deflection and its associated derivatives. The 
latter approach was used in some of the references indicated above and in- 
volves no essential difference in solution techniques. 


Theory 


A. Fundamental Solutions 


1. Simply Supported Plate 
The plate is referred to a cartesian coordinate system (Figure 1). The 

positive z-direction is assumed to be taken into the plane of the paper. 

a) General Loading 


This problem has been completely explored(1) and need not be considered 
here. It is only important to note that the stress, deflection and edge reac- 
tions can be found for this case. 


b) Temperature Gradient 


The problem of a simply supported rectangular plate subjected to a tem- 
perature distribution varying linearly through the thickness has been solved. (1) 


If the plate is referred to the coordinates of Figure 1, the solution takes the 
form 
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Fig. 1. Simply Supported Plate 


t > 0 implies that the deflection will be in the positive z direction. 
Using the deflection function given by equation (1), the distributed line re- 


actions can be determined.(1) This will be done only for VAN " but can be &. 


generalized for each line loading. Thus 


[Max (2% Mey) @) 


cae Differentiation of each infinite series can be performed termwise since each 
successive resulting series converges uniformly(8) in any chosen closed sub- 
interval of the interval o<y<b. Thus 


cosh gm (Y4- 8) (3) 


Cosh an 


mei,3-- 


The expression (3) represents the line reactions at all points except at the i 
corners. The corner reactions can be determined from equilibrium consider- ‘a 

< ations along each edge, and will not be discussed further here. . 

The Fourier sine expansion of the line reactions along the edges will be 

ag obtained for future use. Due to the non-uniform convergence(§) of the series 
representing this reaction V over o<y<b the coefficients of the Fourier sine be 

series cannot be determined by interchanging the order of summation and in- - 
tegration. To by-pass this difficulty, use is made of the fact that the series + 

(3) converges uniformly in any closed subinterval of the region o<y<b. + 
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Therefore the function [%. ] may be approximated to any prescribed de- 


x=a 
gree of accuracy over a given sub-interval by replacing the right hand side of 
equation (3) by a finite series (see Appendix A). All of the ensuing steps for 
the determination of the coefficients of the Fourier series are then valid. 
Thus equation (3) can be replaced by 


cosh Gm ( Y- 8) 


COSH Am 


m= 


This can be represented in the open interval o<y<b in the terms of the sine 
series as, 


where the A, are determined from 


2 , (6) 


Combining equations (5) and (6) and interchanging the order of the finite sum- 
mation and integration there results, 


Aw - (7) 
m by 1+ (any 


2. Edge Loaded Rectangular Plate 
The deflection of a rectangular plate (Fig. 2) loaded along one edge and 
simply supported along the other three edges will now be discussed. 
Subject to the basic differential equation and the boundary conditions of 
zero moment everywhere and zero deflection along the supported edges, and 
in addition 
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Fig. 2. Edge Loaded Rectangular Plate. 


along the loaded edge (x=a), where Ap,’ are the Fourier coefficients asso- 
ciated with the line loadings there results, 


co 


Ws F, ( cosh sinh sin (10) 
ma 1,2, 
From equations (9) and the definition of F,,, the distributed line force along ea 


the edge A (x=a) is 
1 11 
] = Fat sin (11) 


Let v represent the line load along edge j due to a line load along i, where 


positive V implies the loading in the direction of positive z on each edge. 
Then, with this convention, (1) 


x) 
7 

% 
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D[ + (2-V) Wyxx 
= w + (2- Vv) Wyyy] , 
= + o| Wyyy + (2- Vv) Wyarl, 
By direct calculation, from equations (10) and (12) 


Vi = D) En A sin y (13) 


B. General Solution-Elastically Supported Plate 


The general problem of a rectangular plate with four edges elastically 
supported, subjected to an arbitrary load or thermal gradient, will not be in- 
vestigated. This is solved by the successive superposition of the three types 
of structures shown in Figures 3(a), to 3(c) inclusive. Two conditions must 
be met at each edge of the plate. These are: 


1. The deflection of the edge of the plate must be identical with the deflec- 
tion of the supporting beam (condition of compatibility). 


2. The distributed line reaction along each edge must be supplied by the 
beam, so that the corresponding load on the supporting beam must be equal 
and opposite (condition of equilibrium). 


These conditions are satisfied in equations (14) and (15) below, 
Compatibility: 


( Weiate); 
( A, BC, 0) 
Equilibrium: 


+V; ray, +V. (Vecam); z O 


Where all loads are considered positive in the direction of positive z. 
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The subscript i refers to the edge being considered and the superscript 


refers to the loaded edge which defines the structure associated with 
figure 3(b). 


The effects of a line load along edge A can be obtained from equations (10), 
(11) and (13). 

The effect of a line load along B may be obtained most directly from 
Figure 3(b) and the associated equations (11), and (13), if (x, y) and (a, b) are 
replaced respectively by (y, x) and (b, a). 

With reference to Figures 3(b) and 4, the reactions and displacements due 
to a line load along C [ x', y coordinates ] will be identical in form to those 
obtained for the plate loaded along A [ x, y coordinates . The x' coordinate 
may be eliminated through the transformation x' =a - x. 

The effect of a line load along D may be obtained from the result for a line 
load along C in exactly the same manner as B was obtained from A. 


Fig. 4. Line Loaded Plate Referred to x' y Coordinate System. 


This leads to 


w*= > Cosh - sinh sin 


Fe cosh G.y-R,, sinh sing x 


(16) 
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Ve = y cosh g_y} 
= 2(2-v)]sinh g ( (b-y) cosha(o)} 


= sinh (b-y)-(i-v)q(b-y) coshq(b-y)} 


where 


w! = Deflection of plate due to a line load along edge i 


F = Coefficients to be determined 


In determining the Fourier sine expansion for the line reactions given by 
equation (18), it is assumed that the order of summation and integration can 
be interchanged. After simplification, there results 
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As indicated in section A, the downward positive line loads along the edges 
of the arbitrarily loaded simply supported plate can be determined and ex- 
pressed in terms of a Fourier sine series, 


Vase =), Cm SIM 


Veue -) sin q x 


(20) 


As an example, the quantities A,, are determined from equation (7) for the 
simply supported plate subjected to a thermal gradient. The quantities By, 
Cm, Dm, may be determined analogously for this problem. 

From equation (14), considering the problem types being superposed, the 
beam deflection can be written as 


(21) 


(i= A,B,C, D) 


where the quantities on the right hand side represent the it edge deflection 
for a plate simply supported on three sides, and edge loaded on the ith side. 
The values of Wi can be obtained from equations (16) and written as 


m 


sinh ¥,,@ 


For i= A,C 
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x 


The beam deflection and loading are related by, 


From equation (22), the bracketed term in equation (24) becomes 


Expressing (Er) = (EL). \+ gin) ] » and assuming that the 


order of integration and summation can be interchanged, equation (25) can be 
written as the Fourier series 


d? 
et & sin 
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om 


* qm) sin sin dn 


For a beam of constant cross-section, the infinite summation in equation (27), 
is identically zero since g(7) = 0. 
Combining equations (24) and (26) there results 


and for i= B, D 
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= 
Ag 
2 2 4 
4 
am _ : 
# 
(26) 
where 
and 
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Veeam = - Gy Sin (28) 
em 


Beam 


Substituting equations (17), (19), (20) and (28) into equation (15) and replac- 
ing the quantities C 1 by their definitions, the basic equations for the deter- 
n i 
mination of the unknown coefficients F, become, 


my am tn + mit on 
(29b) 
A 2 2 ¢ 
2 20-0 D 


om mn (29d) 


Once these coefficients are determined, the deflection can be evaluated by 
superposition as 


W = Weg t W4A+W HW 


where the wi are obtained from equation (16). 

In the special case of constant (EI)}, the coefficients ct. are zero so that 
the associated summations disappear in equations (29a) to (29d) inclusive. 
C. Summary of Method for Obtaining Solution 


The following steps are suggested for the solution of the elastically sup- 
ported rectangular plate subject to mechanical or thermal loads. 


1. Calculate the reactive force along each of the supports for the simply 
supported plate acted upon by the particular mechanical loading or thermal 
gradient in question. 


2. Expand these functions in a Fourier sine series in order to determine 
Am, Bm, Cm, Dm in equation (20). 


3. Evaluate Bm (1), Bm,n(2), Kni, Cin» Tn» Sn» Wn, and tp. 
B D 
4. Solve equations (29) for nil Fm FC, and Fn: 
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5. Substitute these coefficients into equations (16) and then evaluate w(x, y) 
using equation (30). 


6. Evaluate all required stresses by using the proper combinations of the 
derivatives of the deflection. (1) 


D. Special Cases 


Several of the possible combinations are discussed to illustrate the details 
of the solution technique. 


1. Two Opposite Edges Simply Supported and the Other Edges Elastically Sup- 
ported on Constant Rigidity Beams 

Edges B and D ee ah simply supported by letting the associated beam 
rigidities and hence KB and KD approach infinity. Using these values in 
equation (29b) and (29d) it follows that F,B = F,P = 0. Equations (29a) and 
(29c) then reduce to, 


A 


Solving for FA and FS there results 


Equation (30) then reduces to 


where these separate deflections wh and w© can be obtained from F,4 and 
F,© together with equation (16). 
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2. Two Opposite Edges Simply Supported on Constant Rigidity Beam and the 


Other Two Edges Free 
The equations (31) to (34) inclusive apply for this case where 


« £2, 
so that 
Kn = O 
therefore 
mAmt CoV im | 
where 


a= D[r?- 


(35) 


(36) 


(37) 


(38) 


3. Three Edges Simply Supported and the Fourth Edge Elastically Supported 


on a Beam of Constant Rigidity 


Let the edge A be the elastically supported edge, then KB, KS, and KD 


approach infinity. From equations (29b), (29c) and (29d), there results 


Equation (29a) then yields 


A_ Am 
O(r,+ KS) 
Thus 
W= 


4. Three Edges Simply Supported and the Fourth Edge Free 
For this case Egl, = 0, so that KA = 0. Therefore 


(39) 


(40) 


(41) 


(42) 


| | 
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5. Square Plate, Subjected to a Thermal Gradient, Three Edges Simply Sup- 
ported and the Fourth Edge Free 

This is a special case of 4. The deflection of a simply supported plate is 
given by equation (1). The deflection wA is obtained from equation (16), using 
FA from equation (42), From equation (41) the stresses can be evaluated(1) 
to be 


(43) 


= 


* sinn 


sinh mir 


SIT 


COSH 


An 9..X Cosh Coth mt sinh g.x 
(3+V)Cosh mT, (I-V)_ mir 
sinh mir 


Ma’ 4 mir COSN Ky» 
ms: odd (48) 


het, cosh 9, x+ {> coth mmsinh 4, 


(3+V)COSh MTT + 
sinh mir 


(47) 


The results for deflection and bending moments are plotted nondimension- 
ally in Figures 5 to 7 inclusive. The plate with four sides simply supported 
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is designated by “Simple Support” and that with the one free edge by “Free 
Edge.” Both are included for purposes of comparison. 

As is well known,(1) the value of T along the free edge, as determined from 
shear force Q, need not vanish. This results from the use of the elemen- 
tary theory and the necessary combining of the shear force and twisting mo- 
ment boundary condition. The largest shear stress occurs at the free edge, 


and for the square plate it can be shown, by calculation, to be related to the 
largest bending stress by 


Zh 


Thus even relatively thick plates with a thickness to length ratio of 1/10 
would give a maximum shear stress amounting to, at most, 7% of the direct 
stress. For thinner plates this shear stress is correspondingly less. In ad- 
dition these maximum stresses do not occur at the same location, so that the 
shear stress becomes negligible in design considerations. 


Appendix A: Finite Series Approximation 


Consider a rectangular plate subjected to a loading q(x,y) and reactions as 
shown in Figure A-1, 


Fig. A-1. 


The deflection solution for this general problem can be obtained from the 
solutions for the three superposeable loadings shown in Figures A-2a to 
A-2c. Figure A-2 (a) illustrates the loaded plate with corner reactions only. 
Figure A-2 (b) illustrates the symmetrical component of the line load. 
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Figure A-2 (c) illustrates the anti-symmetrical component of the line load. 
It is necessary to investigate the error introduced by finite series approxi- 
mation to the line load. This error will affect the plate of Figures A-2(b) 
and A-2(c) only. The symmetrical case will be discussed and the anti- 
symmetrical situation follows analogously. 

Consider the distributed symmetrical line load on a representative edge 
together with its equilibrating reactions at the corners, Fig. A-3. Divide the 
load and associated reactions into 2 parts as shown in Fig. A-4, where T = 

: T' + T", and each sub-system is in equilibrium. 

The loads acting over the central region in Figure A-4 (a) can be replaced 
by the statically equivalent system (Fig. A-5). The moment 


M ‘a 
can be made as small as desirable by decreasing 6. Furthermore, from St. 
Venant’s principle, the effect of the self equilibrating force system associated 
with T' is negligible except in the neighborhood of the corners. 

By taking a sufficiently large number of terms, the actual line load may be 
approximated uniformly well in the central sub-region, Figure A-4 (b), and 
the interval of error (6) correspondingly approaches zero. The associated 
Fourier sine series expansion should always include enough terms to repre- 
sent the exact line loading well, except possibly in the subregion of 
Fig. A-4 (a) of width §. This can be always ascertained in advance of obtain- 
ing the complete solution. Therefore the use of the finite series will lead to 
stresses and displacements which agree sufficiently well with that obtained 
for the exact line loading except possibly in the neighborhood of the corners. 
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INCREMENTAL COMPRESSION TEST FOR CEMENT RESEARCH@ 


Corrections by A. Hrennikoff 


CORRECTIONS. —The following illustrations were inadvertently reproduced 
in Proc. Paper 1604 in a size too small for clarity. To facilitate study of 


these figures, they are reproduced herewith. 


Comptession Test July 13, 1956 


Spetimen $8 2'9 x 4° tony Moist 


Age + appr ii months 


Prapertion By wt -24-4, poweotth water cement 


Forwayd Test: Leading) 
Strain 4630 / inch 
Totat | Creep 2155 m-in jin 
Total Elastic Strain main Vin 
Barkward Test) ( Untcading) 
Total: Strain om-in/in. in wonutes 
Tofali Creep | [335 ma/in. 
4 E+ 4150 
Tofal | Elastic Strain 2695 m-m/in. Mean 
Residual Creep not atiowing for 


Spedimen Unbroken 


bey 62 


Note: Paper 1724 is part of the copyrighted Journal of the Engineering Mechanics 
Division, Proceedings of the American Society of Civil Engineers, Vol. 84, No. 


EM 3, July, 1958. 
a. Proc. Paper 1604, April, 1958, by A. Hrennikoff. 
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Compression Test | 


long Dry | (Dried to coast wl 2004) 

Ag e:-Ke most Yor dried For 

ion by wt: 0.003 : BCicement. | 
/ 


Total Strain ( i 


3/39 micro- /inch | | in 7O 
Total Creep 120 m-in /in\ | 


Mean | 4205 “Ys, in| 
Backward Test (Unloading 
Total Sfrain | JOYS infin 


Total 27 musin 
Tatal strain 30/ jntin. | 


Specimen | | 
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i fy 10\, 1957. | 
| Test | 


Specimen 6CIP - 29% 4 ‘Long . Dr (Orved tyconst wt @ 200%) | 
Age - wept moist for 7 months, then dried For. 5 months. 
by wt: 2.003 : water : \cement . 
dest 


ra ld Test (Loadin 


Total Strain ( ing) ince | in 7O min 


Total Evastie m-in7in | 
| Mean E 4/40 in| 
| Backward Test (Unloadin ) | 


Total Creep 30 m- infin 
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